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The  operating  characteristics  of  several  phase-processing  receivers  lor  detection  of  a 
sine  wave  in  Gaussian  noise  arc  deduced  both  analytically  and  by  simulation  and  compared 
with  optimum  processors  over  a  wide  range  of  signal-to-noise  ratios.  Among  the  processors 
considered  are  (I)  the  optimum  phase-processor  for  small  signal-to-noise  ratio  and  known 
sine  wave  frequency,  (2)  the  maximum  likelihood  phase-processor  for  small  signal-to-noise 
ratio  and  unknown  sine  wave  frequency,  (,'j)  the  phase-difference  processors  for  known  or 
unknown  sine  wave  frequency,  and  (4)  a  phase  processor  that  measures  the  scatter  ol  phase 
samples  about  the  best-fitting  straight  line.  It  is  found  that  using  phase  samples  alone  is 
much  better  than  using  amplitude  samples  alone,  and  is  not  much  less  effective  than  using- 
amplitude  and  phase  samples  together  when  the  number  of  samples  is  large  and  the  signal- 
to-noise  ratio  of  each  sample  is  small;  loss  of  amplitude  information  causes  about  J.  dB 
degradation  in  m;my  important  situations. 
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receivers  for  detection  of  a  sine  wave  in  Gaussian  noise  are 
deduced  both  analytically  and  by  simulation  and  compared  with 
optimum  processors  over  a  wide  range  of  signal-to-noise  ratios. 
Among  the  processors  considered  are  (1)  the  optimum  phase- 
processor  for  small  signal-to-noise  ratio  and  known  sine  wave 
frequency,  (2)  the  maximum  likelihood  phase-processor  for 
small  signal-to-noise  ratio  and  unknown  sine  wave  frequency, 
(3)  the  phase-difference  processors  for  known  or  unknown  sine 
wave  frequency,  and  (4)  a  phase  processor  the*  measures  the 
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ber  of  samples  is  large  and  the  signal-to-noise  ratio  of  each 
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1  dB  degradation  in  many  important  situations. 
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DETECTION  CAPABILITIES  OF 
SEVERAL  PHASE- PROCESSING  RECEIVERS 


1.  INTRODUCTION 


Optimum  detection  of  a  finite-duration  sine  wave  of  known  frequency  in 
Gaussian  white  noise  involves  the  use  of  a  matched  filter.  The  matched  filter 
output  is  sampled  at  the  appropriate  time  instant  if  phase-coherent  reception 
is  used,  whereas  the  envelope-detected  output  is  sampled  if  the  received  signal 
phase  is  unknown.  The  sample  is  then  compared  with  a  threshold  to  decide 
whether  a  signal  is  present  or  absent.  Three  major  practical  limitations  that 
preclude  the  use  of  these  processors  arc  that  either  (1)  fhe  duration  or  location 
of  the  sine  wave  may  not  be  known,  (2)  the  medium  frequency-spread  may  not 
be  known  (as  for  example  in  a  fading  medium),  or  (3)  the  bank  of  filters  required 
to  match  unknown  Doppler  shifts  may  be  too  large  or  too  expensive  to  build,  es¬ 
pecially  for  a  long  duration  sine  wave  and  a  multibeam  receiving  system.  An¬ 
other  disadvantage  of  such  processors  is  that  unknown  noise  levels  require  con¬ 
stant  monitoring  and  threshold  adjustment  in  order  to  realize  a  specified  false 
alarm  probability.  (The  losses  associated  with  this  thresh  )ld  adjustment  for  a 
phase-coherent  receiver  are  documented  in  reference  1.  ) 

What  is  desired  is  a  baric  of  fewer  filters  of  broader  bandwidth  than  the 
unknown  matched  filter,  accompanied  by  an  adaptive  processor  that  combines 
the  appropriate  number  of  broadcr-filter  outputs  for  near-optimum  detection 
and  possesses  a  constant  false  alarm  rate  (CFAR)  capability  in  unknown  noise 
levels.  One  possible  way  of  partially  fulfilling  these  requirements  is  based  on 
the  following  observation:  Even  when  the  output  signal-to-noisc  ratio  (SNR)  of 
a  narrowband  filter  has  dropped  to  the  point  where  signal  amplitude  indications 
have  vanished,  a  time  history  of  the  output  phase  samples  continues  to  provide 
straight  line  patterns  that  can  be  discerned  by  a  human  observer.  Since  the  ob¬ 
server  automatically  searches  out  and  averages  the  appropriate  length  of  time 
history  from  the  phase  chart  where  the  frequency  line  is  stable,  he  is  effectively 
performing  a  function  approximating  coherent  processing  over  the  appropriate 
(but  unknow  i  a  priori)  time  interval;  this  function  can  be  interpreted  as  adaptive 
narrowband  filt<  ,  :  g.  The  amount  of  additional  SNR  that  can  be  gained  from  this 
phase-processing  technique  depends  on  the  received  line  stability",  the  filter 
bank  bandwidths,  the  human  observer's  capability',  and  the  operating  point  (i.  e.  , 
false  alarm  and  detection  probabilities). 
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The  advantages  offered  by  this  proposed  phase-processing  technique  could 
be  realized  (and  surpassed)  using  a  filter  bank  of  the  appropriate  width.  However, 
since  the  appropriate  width  is  not  known  a  priori  (and  will  probably  vary  with 
frequency,  location,  environment,  and  local  conditions)  the  proposed  suboptimum 
technique  is  definitely  worthy  of  consideration  and  will  be  investigated  in  detail 
here  to  determine,  quantitatively,  the  amount  of  improvement  achievable  in 
practical  situations.  Actually,  several  phase-processing  techniques  will  be  in¬ 
vestigated  and  compared  with  optimum  processing  in  order  to  determine  the 
limits  of  the  capabilities  of  each  The  CFAR  capability  will  also  be  discussed 
when  appropriate. 

Section  2  presents  a  definition  of  the  problem  to  be  considered  and  the 
forms  of  processing  to  be  compared  for  use  with  both  known  and  unknown  sine 
wave  frequency.  Section  3  provides  quantitative  evaluations,  obtained  both  analyt¬ 
ically  and  by  simulation,  of  the  performance  of  the  various  processors.  Section  4 
contains  conclusions  about  the  relative  merits  of  each  processor.  The  deriva¬ 
tions  of  the  processing  forms  are  presented  in  appendix  A;  useful  geometrical 
interpretations  of  some  of  the  processors,  for  both  known  and  unknown  sine 
wave  frequency,  are  derived  and  discussed  in  appendix  B;  derivations  of  proc¬ 
essor  performance  capabilities  are  collected  in  Appendix  C;  and  the  simulation 
procedure  and  program  is  documented  in  appendix  D.  For  ease  of  cross  refer¬ 
encing  the  appendixes  are  arranged  such  that  particular  processors  have  the 
same  section  numbers  in  each  appendix;  e.  g.  ,  Processor  V  for  known  signal 
frequency  is  treated  in  sections  A.  1.  5,  B.  1.  5,  and  C.  1.  5. 

2.  PrtOBLEM  DEFINITION 
2.  1  KNOWN  SIGNAL  FREQUENCY 

The  received  signal  (if  present)  has  amplitude  PQ,  phase  i PQ,  and  fre¬ 
quency  fQ.  The  sine  wave  frequency,  fQ,  is  assmned  know,  in  this  subsec¬ 
tion;  tliis  restriction  is  removed  in  Section  2.  2.  The  sine  wave  amplitude,  Po, 
is  also  assumed  known.  As  will  be  shown,  however,  none  of  the  processors 
considered  require  or  use  knowledge  of  P(),  i.  e.  ,  they  are  uniformly  most 
powerful  with  respect  to  PQ  (see  reference  2,  pp.  88-92).  Therefore,  this 
restriction  can  also  be  removed.  Various  assumptions  about  the  orient  of  knowl¬ 
edge  of  sine  wave  phase,  ^0,  are  made  below.  The  received  signal,  s(t),  is  ex¬ 
pressed  mathematically  as  a  function  of  time,  t,  as 

s(t)  -  PQ  cos  (2*f^t  +  iG),  t  *  T  ,  (1) 
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where  T  is  the  signal  duration.  This  model  allows  for  multipath  in  the  medium, 
but  not  for  fading.  More  precisely,  any  fading  must  be  slow  enough  so  that  no 
significant  amplitude  or  phase  change  occurs  during  the  time  interval  T 

2 

The  received  noise,  n(t),  has  a  variance  a  and  is  stationary,  narrow- 
oand,  and  Gaussian.  Although  a %  is  assumed  known,  several  of  the  phase 
processors  will  not  require  or  use  this  knowledge.  The  received  waveform, 
r(t),  is  then  given  by 


r(t)  = 


s (t)  t  n(t),  signal  present 

i 

n(t),  signal  absent 


The  processes  described  in  (1)  and  (2)  are  considered  to  be  the  outputs  of 
a  narrowband  filter  (NBF)  of  some  convenient  bandwidth  in  a  receiver  preproc¬ 
essor.  These  NBF  outputs  a»e  sampled  at  time  intervals  approximately  equal 
ro  the  inverse  filter  bandwidth,  so  that  the  in-ohase  and  quadrature  output  noise 
component  samples  are  statistically  independent.  (If  the  received  h,nc  is  not 
centered  in  the  NBF,  P  is  smaller  than  it  would  have  been  for  a  centered  fil¬ 
ter;  thus,  mismatch  of  filter  center  frequency  is  easily  incorporated  in  the  size 
of  r().  Of  course,  for  known  signal  frequency,  this  mismatch  can  be  avoided  ) 

The  received  waveform  (for  signal  present)  can  be  represented  in  complex 
envelope  notation  (see  reference  3,  chapters  I  and  II  or  reference  4,  appendix 
A)  as 

r(t)  =  Re  j  exp  (i2 tt fQ  t)  [PQ  exp  (i  ^Q)  +  n(t)]|  ,  (3) 

where  fQ  is  also  used  for  the  noise  center  frequency,  and  noise  complex  en¬ 
velope  n(i)  is  given  in  terms  of  its  real  in-phase  and  quadrature  components 
according  to 

n(t)  -  nc(t)  +  ins(t)  .  (4) 

A  sample  of  the  complex  envelope  of  the  received  waveform  at  time  tk  is  there¬ 
fore  given  by 

P  exp  (i  vto)  +  n(tk)  -  Pt)  cos  (tQ)  *  r.c  (t,  )  4  i  [PQ  sin  (^Q)  t-  ns(tk)] 

(5) 

=  *k 1 1  yk  =  Rk exR  9  k)  » 
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where  is  the  k-th  in-phase  sample  and  is  the  k-ch  quadrature  sample  of 
the  received  waveform.  Alternately,  Is  the  k-th  amplitude  sample  and  6^ 

is  the  k-th  phase  sample.  For  M  samples  available,  1  £  k  £  M,  the  problem 
we  will  address  is  how  j  R^|  ^  or  j  #}.j  ^  or  both  should  be  processed  for 
near-optimum  detection  capability  when  different  degrees  of  knowledge  of  signal 
phase  i j,  and  frequency  fr  arc  available,  and,  in  some  cases,  for  small  SNR. 

The  complex  sample  x  ^  •  i  yg  can  also  be  considered  to  be  the  output  of 
the  particular  frequency  bin  at  £q  of  a  fast  Fourier  transform  (FFT)  applied 
directly  to  the  received  waveform  prior  to  any  prefiltering.  The  time  between 
samples  is  then  the  time  interval  over  which  the  FFT  is  conducted.  The  follow¬ 
ing  investigation  is  concerned  with  the  way  in  which  the  sequence  of  FFT  outputs 
should  be  combined  and  how  well  the  various  combinations  perform. 

Seven  processors  that  wi  1 !  be  investigated  for  their  detection  capabilities 
are  listed  in  table  1.  'I  he  decision  variable  for  each  is  compared  with  a  thresh¬ 
old.  If  the  th'-eshold  is  exceeded,  a  signal  is  declared  present;  otherwise  it  is 
declared  absent.  (For  Processor  VI,  the  reverse  comparison  is  made.  )  The 
derivations  of  the  processors  in  tabic  1  are  given  in  appendix  A,  along  with  the 
assumptions  necessary  to  make  each  valid  or  optimal. 

Table  1.  Seven  Signal  Processors  for  Known  Signal  Frequency 

Knowledge  Assumed  About  ,  .  ,,  ,  , , 

Processor  „  ,  ,  Samples  Used  Decision  Variable 

Signal  Phase  i0 _ 

Amplitudes!^  ™ 


nmpiibuuco  . 

Phases  f>k  He  |e_i  k  S 


Amplitudes  H^  M 

Constant  but  Unknown  Phases  9k  Rg  e*®k 

k  =  1 


EC, 

R, 

k  1 


Phases  0.  Re 


Constant  blit  Unknown  Phases  0, 


Constant  but  Unkn  >wn  Phases  0, 


£  ei0k) 
V  k=l  ) 


M 

L  ci0k 
k  -  1 


M  v2 


-  z  9l't-  L  0 
Mk  i  \Mk  i  7 


VII  Constant  but  Unknown  Phases  0.  Rc<  Y]  exp(l0.  -  i0L._ i) ; 

k  Ik  2 


3  •'  •*' Jv/Xi 


yvryrTr?}"'  I  y>VT^?WJI^saT!p  vyry.*  yC^T* 


V1 WJW^JJgS?  fSS!®^  rrw-  ggiwjflK 


TR  4529 


i 


I’he  first  processor  m  table  1  assumes  knowledge  of  signal  phase  i/-0  one1  uses 
both  amplitude  and  phase  samples  j  Rj, ! ,  jo^l  .n  its  decision-making  procedure. 

It  is  the  0|)timum  phase-coherent  processor;  no  other  is  as  effective  using  the 
available  samples.  A  geometrical  interpretation  of  its  decision  variable  is  that 
an  M-step  random  walk  in  the  complex  plane  is  taken,  the  k-th  step  is  of  length 
R|.  and  angle  0p,  and  the  resultant  sum  is  rotated  by  the  known  amount,  -  \pQ 
radians,  before  its  real  projection  is  compared  with  a  threshold. 

Processor  11  does  not  assume  knowlr  dge  of  signal  phase  \pQ.  In  this  case  the 
length  ol  tlie  total  random  walk  is  compared  with  a  threshold,  k  is  seen  that  if 
|  9  I .  |  tend  to  cluster  around  a  particular  Value  (4/0),  the  length  of  die  total  walk 
will  be  large  and  .signal  detectability  will  be  enhanced.  It  is  shown  in  appendix  A 
that  Processor  LI  is  optimum  v.hen  the  a  priori  probability  density  function  (PDF) 
of  sigi  il  phase  is  uniform  over  2 n  radians.  Also,  Processor  II  results  when  the 
maximum  likelihood  (ML)  estimate  of  signal  phase  is  used,  instead,  in  a  genera' - 
izcu  likelihood  ratio  (LR)  test.  Thus  two  different  approaches  to  the  treatment  of 
signal  phase  yield  the  same  processor. 

In  case  III  we  assume  that  signal  phase  is  independent  every  sample.  The  opti¬ 
mum  processor  iD  this  case  is  the  inl0  processor  (see  appendix  A  or  reference 
2).  An  approximation  to  this  processor  for  small  SNR  is  afforded  by  the  decision 
variable  in  table  1;  it  discards  all  phase  sampl  information.  An  advantage  of  the 
approximation  is  that  the  decision  variable  is  independent  of  P  ,  whereas  the 
inIQ  processor  requires  knowl edge  of  both  PQ  and  a  for  its  realization.  Thus 
the  approximation  is  uniformly  most  powerful  with  respect  to  P  . 

Both  |  Rk  |  and  |9|  [  were  available  for  use  by  Processors  I,  II,  and  III  for 
every  sample  instant.  Since  increased  noise  levels  cause  an  increase  in  the  value 
of  the  decision  variable  in  these  processors,  it  is  readily  apparent  that  none  of 
them  cam  possibly  be  a  CFAR  receiver.  In  contrast,  the  last  four  processors 
shown  in  the  table  are  lestrictcd  to  operating  on  phase  samples  only,  ;ind  will  pos¬ 
sess  CFAR  capability  for  a  fixed  known  valueof  M,  the  number  of  samples.  For 
example,  Processor  IV-assumes  knowledge  of  signal  phase  \pQ,  and  its  decision 
variable  is  identical  to  that  of  Processor  I,  except  that  the  length  of  each  individual 
random  walk  is  equal  to  unity.  It  is  shown  in  appendix  A  that,  for  small  SNR,  Proc¬ 
essor  IV  is  the  optimum  processor  operating  solely  on  phase  samples;  i.  e.  ,  no 
other  processor  restricted  to  employing  only  |9j.  |  can  do  better  for  small  SNR. 

Processor  V  does  not  require  that  \pQ  be  known,  but  docs  require  that  ^ r 
be  constant  during  all  M  samples.  The  resultant  processor  is  identical  to  Proc¬ 
essor  II  except  that  the  'length  of  each  individual  random  walk  is  again  restricted 
to  unity.  It  is  shown  in  appendix  A  that  Processor  V  is  optimum  when  the  a  priori 
PDF  of  the  signal  phase  is  uniform,  the  SNR  is  small,  and  only  phase  samples  are 
available  for  decision- malting.  Also,  Processor  V  results  when  the  SNR  is  small 
and  the  ML  estimate  of  signal  phase  is  used,  instead,  in  a  generalized  LR  test. 
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Processor  VI  attempts  to  characterize,  mathematically,  what  happens  when  a 
human  observes  a  phase-versus-time  chart  constructed  from  Specifically, 

the  best  horizontal  line  is  fitted  to  the  phase-time  chart  (for  known  signal  frequency 
fy),  such  that  the  average  squared-error  over  M  samples  is  minimized.  The  re¬ 
maining  scatter  about  the  best  straight  line  is  used  as  a  decision  variable  (sec 
appendi  (  A  for  derivations).  The  samples  need  not  be  equi-spar  jd  in  time.  There 
is  an  ambiguity  inherent  in  Processor  VI  because  of  the  nonuniqueness  of  phases 
within  multiples  of  2n ;  this  problem  andmethodsof  circumventing  it  are  discussed 
in  section  3. 

Processor  VII  is  aphase-differenceprocessorthatisnotoptimal  in  any  sense. 

It  is,  however,  a  pracLeal  compromise  that  can  be  resorted  to  if  the  phase  of  the 
signal  is  expected  to  make  an  an.  apt  change  somewhere  in  the  middle  of  the  sequence 
of  M  samples,  or  if  the  phase  drifts  significantly  over  the  total  sequence  length. 
An  abrupt  phase  change  drastically  deteriorates  the  performance  of  all  the  other 
processors  in  table  1  (except  Envelope -Processor  HI)  because  the  coherent  addi¬ 
tion  can  become  destiuctive  rather  than  constructive.  However,  for  Processor  VII, 
an  abrupt  phase  change  merely  causes  several  noisy  samples  in  the  coherent  addi¬ 
tion.  Also,  the  samples  need  not  be  o  yii-spaced  in  time.  Processor  Vfl  has  dis¬ 
advantages  that  will  become  evident  when  its  detection  capability  is  investigated 
in  section  3. 

Idle  derivations  for  Processors  I  through  V  that  arc  presented  in  appendix  A 
assume  that  the  additive  noise  is  Gaussian  and  that  the  samples  are  statistically 
independent.  However,  these  processors  can  be  derived  without  either  of  the  re¬ 
strictive  assumptions,  namely,  by  derivations  that  are  based  on  a  geometric 
interpretation  of  the  envelope  and  phase  samples  as  presented  in  appendix  B. 

Derivations  of  the  detection  and  false  alarm  j  robabilities  for  the  sevenproc- 
essors  are  given  in  appendix  C.  Where  possible,  exact  results,  or  a  method  cap¬ 
able  of  leading  to  an  exact  result,  are  presented.  Ir  oome  cases,  only  the  false 
alarm  probability  can  be  calculated  exactly.  In  other  cases,  approximations  to 
the  detection  and  false  alarm  probabilities  arc  given.  I  or  Processor  VII,  a  more 
general  technique  than  is  indicated  in  table  1  is  also  investigated  (see  (C-86)). 

2.  2  UNKNOWN  SIGNAL  FREQUENCY 

For  unknown  signal  hcnuency,  the  performance  of  all  the  processors 
in  table  1  (except  for  Processor  IIP)  will  be  degraded.  The  received  signal 


*Actually,  as  the  unknown  rrequency  shift  of  the  signal  is  made  iarge  enough 
to  decrease  the  signal  output  of  tne  NBF  in  the  receiver  preprocessor,  P  de¬ 
creases  and  performance  degrades.  However,  our  expressions  foc  detectability 
wall  be  in  terms  of  P  itself,  rather  than  in  terms  of  the  signal  mnplitude  in 
the  medium,  the  NBF  transfer- function,  and  the  precise  frequency  shift. 


(i 
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waveform  is  now 


s(t) 


P  cos  (2  7r (f 
o  '  v  o 


fd)  t  H  +q),  t  t  T  , 


(6) 


where  f(j  is  the  unknown  (Doppler)  sliift.  Although  we  again  assume  that  PQ  is 
known,  it  will  not  be  needed  in  any  of  the  tests  to  be  considered.  Since  it  is  un¬ 


realistic  to  assume  that  signal  phase  \pQ  would  be  known  at  some  time  instant, 


while  the  frequency  f^  is  unknown,  there  are  no  analogs  of  Processors  1  and  IV 


here.  Also,  since  the  performance  of  Processor  HI  is  not  affected*  by  a  frequency 
sliift,  we  need  not  reconsider  it  here. 


The  four  processors  of  interest  here  are  listed  in  table  2  and  their  deriva¬ 
tions  are  presented  in  appendix  A.  If  the  decision  variable  exceeds  a  threshold, 
a  signal  is  declared  present;  otherwise  the  signal  is  declared  absent.  (For 
Processor  VI,  the  reverse  comparison  is  made. ) 

Table  2.  Four  Signal  Processors  for  Unknown  Signal  Frequency 


Processor 


VI 


Knowledge  Assumed  About 


Signal  Phase  <p 


Constant  bat  Unknown 


Constant  but  Unknown 


Constant  but  Unknown 


Simples 

Used 


Amplitudes 
Phases  0k 


Phases 


Phases  $ 


k 


Decision  Variable 


i  M 


mai 

f 


Z  Rk  e>p  (iflg  -i2rftk) 

k  =  1 


M 


max 

f 


k  -  i 


exp ( i 6 


-i2xftk) 


1  M  •>  /  i  1V1 

-  TX  -  L 


M 


M  k  1 

12 


\M  k  1 

$  M  - 


For  Processor  II,  the  complex  random  walk  Rj.  exp(i0j.)  is  "unwound"  by 
all  possible  rates  of  rotating  vectors  (see  table  2),  and  the  largest  distance  walk 
is  compared  with  a  threshold.  It  is  shown  in  appendix  A  that  if  the  a  priori  PDF  of 
signal  phase  is  uniform  ai.d  the  ML  estimate  of  signal  frequency  is  employed,  the 
resultant  test  is  as  given  in  table  2.  Alternatively,  the  appendix  shows  that  if  Ml. 
estimates  are  used  for  both  signal  phase  and  frequency,  the  resultant  generalized 
LR  test  is  identical .  Thus  two  approaches  to  signal  phase  yield  the  same  processor, 
and  assumptions  of  small  SNR  arc  not  required. 


kScc  the  footnote  on  page  G. 
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The  physical  intt  rp rotation  of  the  decision  variable  for  Processor  Vin  table  2 
is  identical  to  that  for  Processor  II  except  that  the  length  of  each  component  of  the 
random  walk  is  forced  to-be  unity.  It  is  shown  m  appendix  A  that  this  processor 
results  f'r  small  SNR  when  either  (1)  the  a  priori  PDF  of  signal  phase  is  uniform 
and  the  Mb  estimate  of  signal  frequency  is  used  or  (2)  ML  estimatesof  both  sig- 
n:u  phase  and  frequency  are  used. 

Processor  VI  results  when  the  best  straight  line  is  fitted  to  equi-spaced phase 
samples,  such  that  the  average  squared-error  over  M  samples  is  minimized.  It  is 
felt  that  this  test  approximates  the  way  a  human  observer  somehow  pieces  together 
the  "barber  pole"  effect  caused  by  an  unknown  frequency  when  confronted  with  a 
phase  chart  history  limited  to  the  (-  tt,  7r)  range.  Inpractice,  there  is  a  real  prob¬ 
lem  when  attemptin''  to  apply  the  decision  variable  in  table  2  to  automatic  data 
processing,  because  a  computed  phase  sample  0j.  will  lie  in  the  (-ny  ir)  range  un¬ 
less  a  special  algorithm  is  used  to  "unwind"  the  phase  samples  along  a  straight 
line.  Also,  unless  there  is  human  intervention,  all  algorithms  will  occassionally 
yield  permanent  spurious  j_2ir  jumps  in  phase  that  give  an  unnecessarily  large  value 
to  the  decision  variable  (scatter)  and  can  lead  to  erroneous  signal -absent  statements. 

If  the  \0\i\ could  be  unwound  properly,  the  performance  of  Processor  VI  for  un¬ 
known  signal  frequency  would  be  identical  to  that  of  Processor  VI  for  known  signal 
frequency  and,  therefore,  independent*  of  the  exact,  frequency  shift,  f(j.  Accordingly, 
no  detailed  operating  characteristics  will  be  given  for  this  processor  in  section  3. 

Processor  VII  is  a  phase-difference  processor.  For  equi-spaced  samples  and 
unknown  signal  frequency,  the  complex  vector  exp  (ifg  -  i0g_^)  tends  to  cluster 
around  an  unknown  phase  for  large  SNR.  The  length  cf  the  total  walk  is  used  as  a 
decision  variable,  since  the  direction  of  walk  is  not  known  a  priori.  The  perform¬ 
ance  of  Processor  Vil  is  independent*  of  the  exact  frequency  sluit  fcj,  which  is  can¬ 
celed  out  of  the  decision  variable  (see  (A-3G)  through  (A— 40)).  However,  the  per¬ 
formance  is  poorer  than  that  for  the  comparable  processor  for  known  signal  fre¬ 
quency.  'file  performance  capability  of  Processor  VII  (see  table  2)  is  given  in 
section  3. 

The  derivations  in  appendix  A  for  Processors  II  and  V  (see  table  2)  assume 
that  the  additive  noise  is  Gaussian  and  that  the  samples  are  statistically  indepen¬ 
dent.  The  resulting  processors  are  derived  in  appendix  B  without  cither  of  these 
restrictive  assumptions  via  a  geometrical  interpretation. 

According  to  the  above  discussion,  performance  results  for  Processors  II, 

V,  and  VTI  in  table  2  will  be  presented  in  section  3.  Analytical  derivations  of 
detection  and  false  alarm  probabilities  for  Processors  II  and  V  arc  too  difficult 
and  have  not  been  attempted.  Analytical  results  fer  the  performance  capability 
of  Processor  VII  are  presented  in  appendix  C,  where,  in  fact,  a  processor  of 
a  more  general  form  is  investigated  (see  (C  —  1. 33)). 


*Sec  the  footnote  on  page  G. 
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3.  RESULTS 


3.  1  KNOWN  SIGNAL  FREQUENCY 


Receiver  operating  cha’-acteristics  (ROCs)  for  the  processors  discussed  in 
the  previous  section  were  determined  by  computer  simulation  using  10,000  inde¬ 
pendent  trials.  Although  the  ROCs  for  Processors  I,  II,  and  III  are  available  in 
the  literature  (e.  g.  ,  in  references  3  and  5),  they  are  included  in  this  simulation 
for  purposes  of  determining  the  credibility  and  accuracy  of  the  simu’  .tion  ap¬ 
proach.  Detnj  s  of  the  simulation  used  are  presented  in  appendix  D. 


In  the  following  discussion,  the  number  of  samples,  M,  is  fixed  at  25  and 
comparisons  between  processors  are  made  on  this  basis.  For  smaller  values 
of  M  of  interest,  the  program  in  appendix  D  can  be  suitably  modified.  For 
larger  values  of  M,  the  computer  time  and  storage  increase  and  the  analytical 
derivations  in  appendix  C  can  be  used,  instead;  in  particular,  the  approximations 
in  appendix  C  improve  for  large  M. 


In  figure  1A,  the  ROC  for  Processor  I  is  presented  for  values  of 

P 


m  ~  — 


ranging  from  0  to  .  8.  The  SNR  corresponding  to  (7)  is,  in  decibels, 


Po/2\ 


10  log  l — J  -  20  log  (m  >  -  3 


This  is  the  SNR  it  the  preprocessor  output,  prior  to  recombination  according 
to  table  1.  In  the  following,  Pp  is  the  false  alarm  probability  and  Pp  is  tie.’ 
detection  probability. 


The  (0,  .1)  range  for  false  alarm  probability,  Pp ,  is  expanded  in 
figure  IB  to  clearly  show  the  ROC  for  small  Pp  in  the  neighborhood  of  .01. 

As  a  check  of  the  accuracy  of  the  simulation  approach,  several  exact  values  of 
performance  have  been  added,  as  given  by  (C-5)  in  appendix  C,  and  arc  denoted 
by  Xs  in  figure  1.  T  my  arc  seen  to  virtually  coincide  with  the  simulated  results. 


The  ROC  for  Processor  II  is  presented  in  figure  2.  Again,  several  exact 
values  (from  (C-9))  have  been  entered  as  Xs,  and  the  agreement  is  very  good 
over  a  wide  range  of  values  of  Pp.  In  comparison  with  Processor  I,  the  de¬ 
gree  of  degradation  suffered  by  Processor  II,  which  lacks  knowledge  of  signal 
phase  \t0,  depends  on  the  exact  operating  point  on  the  ROC.  For  example,  the 
point  (Pp,  Pp)  =  (.  01,  .  2)  requires  m  =  .  3  for  Processor  I,  but  m  =•  .  4  for 
Processor  II,  which  is  a  20  log  (.  4/.  5)  -  2  5  dB  degradation.  Processor  II 


w 


M 
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suffers  a  certain  amount  of  small-signal  suppression;  i.  e.  ,  dete.  ion  of  weak 
signals  (m  <<  11  is  significantly  impaired  in  comparison  with  Processor  I.  This 
is  characteristic  of  all  processors  that  do  not  have  knowledge  of  signal  phase 

'/'o' 

l’he  ROC  for  Processor  III  is  given  in  figure  3.  It  is  seen  that  it  suffers 
even  greater  small-signal  suppression;  e.g. ,  a  value  of  m  =  .  78  is  required 
to  realize  he  (.  01,  .  2)  operating  point,  which  is  8.3  dB  poorer  than  Processor  I. 
As  a  comparison  witn  figure  3,  the  ROC  for  the  optimum  ini  processor, 
(A-15),  was  also  computed  for  the  same  set  of  pa ’ameters;  the  results  were 
imperceptibly  different  from  figure  3  and,  therefore,  are  not  presented.  Ex¬ 
tensive  numerical  results  for  the  performance  of  Processor  III  are  available 
in  reference  5. 

The  ROC  for  the  first  of  the  phase  processors,  IV,  is  given  in  figure  4.  Its 
performance  is  slightly  poorer  than  Processor  I,  which  used  both  amplitude  and 
phase  samples.  For  example  the  (.1,  .9)  operating  point  is  realized  with  m  =  .  52 
for  Processor  I,  but  requires  m  =  .  58  for  Processor  IV;  this  is  a  difference  of 
only  .  95  dB.  It  is  shown  in  appendix  C  (equations  (C-40)  and  (C-41)  et  seq.  )  that 
for  large  M  and  small  m.  Processor  IV  is  10  log  (4/w  )  =  1.  05  dB  poorer  than 
Processor  I.  Thus  loss  of  amplitude  information  does  not  cause  a  great  loss  in 
detection  capability;  rather,  the  information  in  the  phase  samples  is  the  domi¬ 
nant  contributing  factor  to  signal  detection.  Since  Processor  IV  uses  knowledge 
of  signal  phase  </- 0 ,  there  is  no  small-signal  suppression;  this  applies  even 
though  |  R|.  |  have  been  discarded. 

Processor  V  does  not  assume  knowledge  of  signal  phase  ^0,  whereas 
Processor  IV  does.  The  ROC  for  V  is  depicted  in  figure  5.  When  figures  2  and 
5  are  compared,  it  is  apparent  that  Processor  V's  perfonriance  capability  is 
very  similar  to  that  of  II,  which  had  amplitude  information  available  in  addition 
to  the  phase  information.  For  example,  operating  point  (.  1,  .  9)  requires  m  =  .  G5 
for  Processor  U  and  m  =  .  74  for  Processor  V  —  a  difference  of  1.  1  dB.  Equa¬ 
tions  (C-68)  and  (C-69)  et  seq.  show  that  for  large  M  and  small  m,  Processor  V 
is  1.  05  dB  poorer  than  Processor  II.  Thus,  once  again,  loss  of  amplitude  infor¬ 
mation  causes  a  loss  of  1.  05  dB  in  detectability  in  one  limiting  case.  Notice  that 
Processors  I,  II,  IV,  and  V  ill  1  use  the  phase  samples  in  the  same  basic  way, 
i.  e.  ,  according  to  the  addition  of  complex  vectors  exp  (iflj.). 

Processor  VI  makes  an  entirely  different  use  of  the  phase  samples,  and, 
because  it  does,  the  inherent  ambiguity  in  phase  causes  some  problems  in  auto¬ 
matic  processing.  To  illustrate  the  problem,  suppose  the  true  signal  phase 
o  =  0  and  are  defined  in  the  (-*■,  ir)  band.  Then,  for  the  decision  variable 
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iii  table  1,  the  corresponding  ROC  for  in  =  0  (.  1)  .  4  is  presented  in  figure  6A. 

This  ROC  is  poorer*  ;han  that  in  figure  4A;  this  result  agrees  with  the  fact  that 
Processor  IV  is  the  optimum  processor  for  small  in,  phase  samples  alone,  and 
known  signal  please.  On  the  other  hand,  the  ROC  in  figure  6A  is  noticeably  better 
than  that  in  figure  5A  because  Processor  \T  assumed  that  -  0.  whereas 
Processor  V  does  not  use  signal  phase  knowledge  (although  it  does  use  |0]c| 
optimally  for  small  m,  otherwise).  Notice  the  absence  of  small-signal  suppression 
in  figure  GA,  since  i/-n  is  known  and  used. 


In  an  attempt  to  determine  the  effect  of  different  signal  phases,  i/-0,  on  the 
performance  of  Processor  VI,  the  decision  rule  in  table  1  was  applied  in  the 
case  where  i p0  =  *'/2  and  |#ic|  were,  again,  in  the  (-rr,  w)  band  The  resulting 
ROC  in  figure  GB  suffers  notable  degradation  because  |0j.|  with  values  in  the 
(- ir,  -  ir/2)  band  heavily  bias  the  average  phase  estimate,  (A-33),  away  from 
the  true  value  \pQ  =  * /2.  Then  the  scatter  (see  table  1  or  (A— 34))  about  the 
average  phase  estimate  is  large  and  decisions  made  about  signal  absence  are 
frequently  erroneous.  Yet,  near  -ir  should  not  be  allowed  to  contribute 

such  a  large  effect  since,  in  reality,  each  value  0  p  •  2*  is  quite  near  the  true 
value  tt/2;  i.  e.  ,  the  (-  jt/2,  3  x/2)  band  would  be  the  ideal  interval  in  which  to 
define  each  0g  when  \p0  -  *  /2.  But  since  the  true  signal  phase  will  not  be 
known  in  practice,  the  band  to  choose  a  priori  is  unknown.  The  problem  becomes 
magnified  when  \)/0  is  allowed  to  sweep  out  a  full  2-n  >'ange  in  the  10,000  trials 
and  band  (-tt,  n)  is  retained  for  {fljJ.  The  resultant  ROC  in  figure  GC  indicates 
that,  for  certain  ranges  of  Pp,  the  processor  is  poorer  than  a  random  choice. 


A  possible  solution  to  the  automatic  processing  problem  here  is  to  define 
each  0^,  simultaneously  in  several  2?r  bands,  and  compute  the  scatter  that 
has  a  minimum  value  (in  all  bands  considered)  as  the  decision  variable.  Ideally, 
there  should  be  a  continuum  of  bands  (-  n  t  v,  w  i  v),  where  v  ranges  contin¬ 
uously  over  a  2  ir  interval,  but,  since  this  is  impractical,  a  few  well-chosen 
values  for  the  bands  should  be  considered.  In  figure  GO,  bands  (- it ,  n)  and  (0,  2 ir) 
were  both  used  as  swept  out  a  full  2*  range  in  the  10,000  trials.  The  im¬ 
provement  afforded  by  the  addition  of  band  (0,  2ir)  is  marked  (compare  figures 
61)  and  GC).  In  fact,  for  small  m,  the  ROC  in  figui-e  GD  is  almost  as  good  as  that 
in  figure  5A  for  the  optimum  phase  processor  without  knowledge  of  signal  phase 


♦In  (C-83)  et  seq.  it  is  shown  that  Processor  VI  is  .34  dB  poorer  than 
Processor  IV  for  large  M  and  small  SNR  m.  Phis  figure  of  .  34  dB  need  not 
hold  true  over  the  entire  ROC,  however;  in  fact,  since  Processor  IV  is  optimum 
only  for  small  m,  Processor  VI  with  known  could  conceivably  outperform 
Processor  IV  for  larger  m. 
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PQ.  Probably,  use  of  the  three  bands  (-5  jt/3,  jt/3),  (-7r,7r),  and(-Jr/3,  5  tt/3), 
with  cen'  sat  -2  7T./3,  0,  and  2^/3,  respectively,  would  produce  a  ROC  that 
is  very  close  to  that  in  figure  5A;  however,  this  is  not  pursued  here. 

Another  possible  solution  attempted  was  to  "unwind"  the  successive  phase 
samples  09,  ,  9  such  that  J  0  k  ~  ^  k  - 1 1  l71" ;  i-  »  the  phase  sam¬ 

ples  were  allowed  to  drift  to  whatever  value  that  was  necessary,  provided  chat 
adjacent  phase  changes  were  never  more  than  vr  radians.  Then,  ti;°  best  tilted 
straight  line  was  fitted  through  the  data  and  the  remaining  scatter  was  used  as  a 
decision  variable  (see  (A-02)  et  scq.  );  the  resulting  ROC  in  figure  (IE  is  totally 
inadequate. 

The  above  discussion  has  dealt  with  an  attempt  to  employ  automatic  data 
processing  for  variations  of  Processor  VI.  If,  instead,  a  human  observes  a 
phase  chart  of  0^,  02»  •••  ,  versus  time,  limited,  perhaps,  to  the  (-*■,*) 
band,  he  will  probably  piece  together  the  "barber  pole"  effect  or  notice  clusters 
of  phase  samples  about  some  constant  nonzero  value,  and  will  be  able  to  make 
accurate  decisions  about  signal  presence,  even  when  \pQ  is  near  +  tt.  (Two  clus¬ 
ters  would  be  observed  for  \ p0  near  ctt  or  for  \pc  near  -ir.)  It  is  not  known 
how  close  the  human  can  come  to  the  optimum  performance  indicated  in  figure 
5 A;  however,  the  similarity  between  figures  5A  and  6D  for  a  very  simple  deci¬ 
sion  rule  and  two  bands,  leads  to  the  conjecture  that  human  pcrfonnance  should 
be  fairly  close  to  that  shown  in  figure  5A. 


A  phase-difference  processor  for  known  signal  frequency,  which  is  a  more 
general  form  of  Processor  VII  than  is  shown  in  table  2,  will  be  considered  here. 
Specifically,  the  decision  variable  for  this  processor  is  a  weighted  sum  of  phase- 
difference  vectors: 


The  reason  for  the  weighting  power  p  is  that,  as  has  been  observed  above,  am¬ 
plitude  weighting  does  improve  performance  somewhat  (compare  I  versus  IV  and 
II  versus  V).  Both  exact  and  approximate  detailed  analyses  of  (9)  in  terms  of 
Pg>  and  Pp  arc  presented  in  appendix  C  for  general  p. 

The  ROCs  for  p  0,  .  5,  and  1,  arc  presented  in  figures  IA,  7B,  and  7C, 
respectively.  These  figures  show  improved  performance  as  p  increases  from 
0  to  1.  In  fact,  in  appendix  C  ((C- 131)  et  seq. )  it  is  shown  that  the  best  value 
of  p,  for  large  M  and  small  m,  is  unity.  Also,  the  loss  in  detectability  at  p  ~  0 
versus  p  1  for  Processor  VTI  is  1.  05  dB  under  these  conditions;  once  again. 


1.2 


two  processors  that  differ  only  in  the  way  they  use  amplitude  infoi'mation  differ 
by  1.  05  dB  in  performance  for  large  M  and  small  m. 


™r?- ' 


It  must  also  be  noted  in  figure  7  that  the  increment  in  m  (  P0/<r)  is  •  2,  not 

.  1  as  in  all  of  the  previous  figures.  Thus,  the  phase-diffei’encc  processor  (lor 
any  weighting  m)  is  obviously  poorer  in  performance  than  the  optimum  phase 
processor,  V,  which  is  without  knowledge  ci  signal  phase  In  fact,  (C-131) 

shows  the  familiar  nm-dependence  of  Pj-p  thus,  the  phase-difference  processor 
suffers  a  small-signal  suppression  effect.  For  larger  m,  a  direct  comparison  of 
ItOCs  is  necessary.  For  example,  although  the  operating  point  (.  1,  .  78)  requires 
m  1  for  Processor  VII  w  ith  n  0,  Processor  V  requires  only  m  .  (13  —  a 
difference  of  4  dB.  On  the  other  hand,  as  was  noted  in  the  Introduction,  the  phase- 
difference  processor  can  tolerate  an  abrupt  signal-phase  jump  without  disastrous 
effects,  whereas  the  other  processors  would  all  be  very  adversely  affected. 


Another  u'ay  of  presenting  the  operating  characteristics  of  a  processor, 
which  can  be  particularly  useful  for  comparison  of  systems  for  small  values  of 
Pp,  is  to  plot  Pp  versus  m  on  probability  paper,  with  Pp  as  a  parameter 
(e.  g.  ,  sec  reference  3,  figures  IV.  1  and  V,  2).  We  shall  call  these  plots  detec¬ 
tion  characteristics  (DC),  in  order  to  distinguish  them  from  ROCs,  and  present 
them  for  several  of  the  processors,  as  determined  by  both  simulation  and  ana¬ 
lytical  methods.  This  approach  will  afford  corroboration  of  the  analytical  results, 
as  well  as  an  indication  of  the  adequacy  and  accuracy  of  large-M  approximations. 
It  will  then  be  possible  to  use  the  analytical  results  for  other  values  of  M  with 
assurance  of  their  accuracy  and  applicability. 


The  DCs  for  Pi’ocessors  I,  II,  and  III  are  available  in  references  3  and  5. 

The  DC  for  Processor  IV  is  presented  in  figure  8,  where  the  solid  curve  repre¬ 
sents  the  approximate  analytical  result  in  (C-40),  with  exact  threshold  values  as 
given  in  table  C-l.  The  superposed  Xs  denote  the  results  of  the  simulation  de¬ 
scribed  above.  The  analytical  result  is  expected  to  be  in  greater  error*  for  small 
m  than  for  large  m.  The  agreement  with  the  simulation  results  is  quite  good; 
the  standard  deviations  of  the  results  arc  considered  in  appendix  D.  For  small 
Pp,  approximation  (C-4  0)  is  not  as  acc  urate,  especially  for  small  in;  in  fact, 
the  limiting  value  Pjr  10"^  is  not  realized  at  m  ;  0  by  this  approximation, 
even  though  exact  thresholds  were  used.  The  simulation  results  for  P v  -  10  1 
are  not  shown  on  figure  9  because  10,000  trials  do  not  constitute  an  adequate 
base  for  estimation. 

♦The  reason  for  this  behavior  is  that  the  individual  random  variables  in 
(C-19)  are  more  nearly  Gaussian  for  large  m  than  for  small  m.  In  fact,  as 
m  —  0  in  (C-21),  the  PDF  of  s  has  cusps  at  j_  1. 
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The  DC  for  Processor  V  is  given  in  figure  !),  where  the  analytical  result, 
(C-(>3),  is  drawn  as  a  solid  curve;  exact  threshold  values  given  by  table  C-2 
were  used,  rhe  superposed  Xs  again  denote  simulation  results  and  it  can  be 
seen  that  disagreement  is  greates,'  for  small  I  *  |.*  and  small  nr,  here  the  approx¬ 
imation  gives  an  overestimate  of  aUnlnuhlc  performance.  The  characteristic 
small-signal  suppression  ot  process!  rs  operating  without  knowledge  of  signal 
phase  \f/Q,  as  exhibited  hv  the  horl/.o!  la!  .slope  at  m  (),  is  apparent. 

A  detection  characteristic  is  not  presented  lor  Processor  VI  because,  as 
described  above,  multiple  bands  were  not  attempted  lor  the  definition  of  |0|.J. 
However,  it  is  anticipated  that ,  lor  the  reasons  discussed  above,  the  DC  in 
figure  9  would  i)e  a  good  appro\i...al  Pai  to  Dial  lor  Processor  VI. 

Three  DCs  for  Processor  Vil  loi  weighling  muslimi  )t  0,  .5,  and  1,  re¬ 
spectively,  are  presented  In  ligun  lo  <"<  '■  <  > M ■  I'he  Ihresiiolds  for  a  0  and 
H  1,  from  tables  C-.’i  and  (  I,  ara  exact,  and  Hit  thresholds  lor  a*  .5  are 
determined  from  ( C - 127).  The  approximate  I’d,  as  determined  from  (C-129), 
agrees  well  with  the  simulation  results  tor  all  three  eases  of  weighting  constant 
A<  considered,  except  in  figure  IOC  lor  n  1  and  small  Pj,..  The  analytical 
results  in  figure  IOC  for  Pj,.  10“''*  and  10  1  are  pessimistic  estimates  of 
performance  for  small  m,  but  are  apparently  optimistic  estimate's  for  larger  m. 
The  reason  for  the  poorer  quality  of  the  approximation  to  P^  for  small  Pp  in 
figure  IOC  is  that  the  summation  random  variables  in  (!))  are  far  from  Gaussian 
for  g  1  and  large  amplitudes.  On  the  other  hand,  the  extremely  good  agree¬ 
ment  indicated  in  figure  10 B  is  fortuitous;  of  course,  the  reason  the  approxima¬ 
tion  (C-129)  is  perfect  at  m  0  is  that  approximate  thresholds,  as  determined 
by  (C-129)  itself,  have  been  used. 

3.  2  UNKNOWN  SIGNAL  FREQUENCY 

The  ROCs  for  Processors  II,  V,  and  VII  in  table  2  will  be  given,  in  keeping 
with  the  discussion  in  section  2.  Again,  10,000  independent  trials  are  taken  and 
M  =  25. 

Suppose  the  complex  samples  | Rjc  exp  (i 0jt) |  are  taken  every  A  seconds; 
i.  e.  ,  suppose  tg  k  A  in  tabic  2.  This  time  interval  corresponds  to  an  NBF 
of  width  l/A  Hz  in  the  receiver  preprocessor  (or  to  the  time  between  FFT  out¬ 
puts),  as  described  in  section  2.  Thus  the  complete  preprocesso  r  might  consist 
of  a  bank  of  NBFs  (or  a  subset  of  the  FFT  outputs)  extending  over  some  selected 
frequency  band.  For  a  collection  of  M  samples,  the  total  observation  interval  is 
MA  seconds;  the  fundamental  frequency  resolution  possible  o'er  this  total  time 
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interval  is  (MA)  ^  Hz.  Therefore,  if  we  expect  the  decision  rules  in  table  2  tc 
yield  near-optimum  results,  the  search  in  f  must  be  conducted  to  within  at  least 
(2MA)'1  Hz. 

Let  the  actual  frequency  search  increment  be  (N  A )'■*■,  where  N  _>  2M  .  Then 
for  f  fn  E  n/(NA),  the  decision  rule  for  Processor  II  in  table  2  takes  the  ap¬ 
proximate  form 

M 

(10) 


max 

n 


Yj  Rk  exP  (iflk)  exP  ("i  2wkn/N)  I  , 
*  1  | 


where  the  search  over  n  covers  only  the  range  of  expected  frequencies.  The 
two  particular  cases  we  consider  are  (1)  a  "track"  situation  where  the  signal 
frequency  is  known  almost  exactly  and  (2)  a  "search"  situation  where  the  fre¬ 
quency  uncertainty  is  approximately  1 M  Hz.  (For  greater  uncertainties,  an- 
o  ner  NBF  in  the  preprocessor  bank,  or  a  different  frequency  bin  in  the  FFT 
output,  would  pick  up  the  signal. )  Mathematically,  these  two  eases  correspond 
to  |  n  |  <  1  and  j  n  |  <.  N/2,  respectively,  in  (10). 


The  summation  in  (10)  is  immediately  recognized  as  an  N-point  FFT  of  the 
M  nonzero  complex  samples  {lip  exp  (iflg)}.  Thus,  (10)  dictates  a  comparison  of 
the  maximum  magnitude  of  the  FFT  (of  the  time  sequence  j  Rg  exp  (ifl  g)  |  ^)  with 
a  threshold  for  decisions  about  signal  presence.  A-,  N  becomes  very  large,  (10) 
approaches  the  continuous  rule  quoted  in  table  2.  However,  computation  of  (10) 
is  unnecessarily  time-consuming  if  N  is  selected  too  large;  also,  performance 
capability  saturates  for  N  ■  2M.  We  will  consider  N  64  here,  because  of  the 
speed  of  the  FFT  for  powers  of  2. 


The  ItOC  for  Processor  II  in  the  track  situation  is  depicted  in  figure  II;  it 
is  slightly  poorer  than  in  figure  2,  which  shows  the  corresponding  ROC  for  known 
signal  frequency.  The  ROC  lor  Processor  II  in  the  search  mode  is  illustrated  in 
figure  12.  At  operating  point  (.  01,  .  62),  for  example,  approximately  1.  8  dB 
more  SNR  is  required  in  the  search  mode  than  in  the  track  mode.  The  difference 
is  much  greater  for  smaller  m.  where  the  lack  of  knowledge  of  signal  frequency 
causes  a  significant  degradation  in  performance. 

For  Processor  V  in  table  2,  the  only  change  in  (10)  is  to  replace  R.,  by 
unity;  then  all  of  the  above  comments  apply  directly.  The  resultant  ROC  is  dis¬ 
played  in  figure  15  for  the  track  mode  and  in  figure  14  for  the  search  mode. 
Comparison  of  figure  13  (unknown  signal  frequency)  with  figure  5  (known  signal 
frequency)  reveals  a  slight  degradation.  The  degrading  effect  of  the  wider  search 
procedure  is  evident  in  figure  I  I. 
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A  phase-difference  processor  for  unknown  signal  frequency,  which  is  a 
more  general  form  of  Processor  VTI  than  is  shown  in  table  2,  will  be  consid¬ 
ered  here.  Specifically,  the  decision  variable  is  a  weighted  sum  of  phase-dif¬ 
ference  vectors: 


M 


£  K  V-l  exp(i0k-i0k  l) 


(11) 


Both  exact  and  approximate  detailed  analyses  of  ^  1 1 )  in  terms  of  Pjy  and  P,r 
are  presented  in  appendix  C  for  general  n. 


The  ROCs  for  a  0,  .  5,  and  Rare  presented  in  figures  15A,  15B,  and 
15C,  respectively.  These  figures  show  improved  performance  as  m  increases 
from  0  to  1.  Equation  (C- 177)  and  table  C-5  in  appendix  C  show  that  the  best 
value  of  for  large  M  and  small  m  ,  is  unity.  Also,  the  loss  in  detectability 
m  v  0  versus  a  1  is  1.  05  dB  under  these  conditions;  i.  e.  ,  when  amplitude 
information  is  discarded  in  this  proc”  soor,  1.  05  dB  is  lost  for  large  M  and  small 
m.  The  ROC  for  a  1.  5  was  also  computed  and  was  found  to  be  virtually  iden¬ 
tical  to  that  for  a  i;  however,  the  ROC  for  ii  2  was  slightly  poorer  than  the 
ROC  for  ;i  1.  (Neither  of  these  eases  is  presented  in  this  report.  ) 


It  Is  important  to  notice  that  the  increment  in  it.  shown  in  figure  15  is  .  2, 
not  .  1.  Thus  the  performance  of  the  phase-difference  processor,  VII,  for  un¬ 
known  signal  frequency  is  extremely  poor  until  m  approaches  unity';  i.  e.  ,  the 
SNR  at  the  preprocessor  output  must  be  approximately  -3  dB  before  adequate 
performance  obtains  (see  ( 8 ) ) .  (The  degradation  of  Processor  VTI  for  known  sig¬ 
nal  frequency  (see  figure  7)  was  not  as  marked  for  small  m.  )  In  (C-177),  the 
approximate  detection  probability  is  shown  to  have  an  irm-dependence  rather 
than  the  m-dependencc  of  Processors  II  and  V.  Thus,  in  addition  to  the  small- 
signal  suppression  effect  caused  by  the  lack  of  knowledge  of  signal  phase  \jy 0 , 
the  p’nase-diffei  ence  processor  for  unknown  signal  frequency  suffers  further 
degradation  for  small  in  for  all  values  of  weighting  constant  m  in  (11). 

Since  we  were  unable  to  derive  analytic  results  for  the  detection  probabili  ;ies 
oi  Processors  II  and  V  for  unknown  signal  frequency,  no  DCs  are  presented  for 
these  two  cases.  However,  we  are  able  to  analyze  Processor  VR  for  unknown 
frequency;  the  resultant  DCs  for  n  0,  .  5,  and  1  are  presented  in  figures  IB  A, 
1GB,  and  10C,  respectively.  'ITic  thresholds  used  for  u  0  and  n  1  are  exact 
and  are  given  in  tables  C-G  and  C-7,  respectively.  The  thresholds  used  for 
ii  .  5  were  determined  from  approximation  (C-173). 
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The  pronounced  flatness  of  the  DCs  for  m  near  zero  is  expected  in  accord¬ 
ance  with  the  suppression  effect  discussed  above.  Thus,  for  example,  even  at 
m  .  5,  the  detection  probability  in  figure  16A  has  only  reached  .  15  for  the  curve 
labeled  Pp  =  10-*.  This  behavior  is  exhibited  independently  of  the  value  of  ^ . 

The  poorer  quality  of  the  approximation  to  Pp,  in  figure  1GC  for  small  Pp 
is  the  result  of  the  fact  that  the  summation  random  variables  in  (11)  are  far  from 
Gaussian  for  p  1  and  large  amplitudes.  The  analytic  results  are  pessimistic 
estimates  of  performance  for  small  m,  but  optimistic  for  larger  m. 

4.  DISCUSSION 

The  failure  to  use  amplitude  information  to  detect  a  sine  wave  in  noise  causes 
a  degradation  in  performance  of  about  1  dB  for  small  SNR  and  large  M  if  the  phase 
information  is  properly  used.  This  conclusion  is  based  upon  a  comparison  of  the 
following  four  pairs  of  processors: 

a.  Processor  I,  known  signal  frequency,  versus  Processor  IV,  known  signal 
frequency; 

1).  Processor  II,  known  signal  frequency,  versus  Processor  V,  known  signal 
frequency; 

c.  Processor  VTJ,  known  signal  frequency,  a  1  versus  a  0;  and 

d.  Processor  VII,  unknown  signal  frequency,  u  1  versus  n  0. 


On  the  other  hand,  discarding  phase  information  can  lead  to  severe  degradation, 
as  the  ROC  for  Processor  III  attests. 


All  processors  that  lack  knowledge  of  signal  phase  suffer  small-signal  sup¬ 
pression,  regardless  of  their  use  ef  amplitude  and  phase  samples.  However,  for 
larger  SNR,  such  that  useful  detection  probabilities  result,  the  degradation  is 
often  insignificant.  The  suppression  effect  is  more  pronounced  for  the  phase- 
difference  processors  and  is,  in  fact,  severe  for  unknown  signal  frequency. 

The  phase  processors  have  CFAR  capability;  i.  e.  ,  a  threshold  can  be  fixed 
to  realize  a  prescribed  Pp  for  a  given  value  r  M,  regardless  of  the  absolute 
noise  level.  The  threshold  for  many  of  the  pr  essors  can  be  determined  by 
means  of  the  exact  formul  is  in  appendix  C. 

Approximations  to  system  detection  capability  derived  in  appendix  C  can  be 
fruitfully  employed  to  evaluate  processor  performance  for  other  values  of  M  and 
SNR  However,  for  small  values  of  ,M  or  'cry  small  values  of  Pp,  the  approxi- 
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mations  are  less  accurate.  A  guide  as  to  when  approximations  are  valid  is  offered 
by  the  DCs  presented  in  this  report.  When  simulation  becomes  expensive,  in 
terms  of  time  and  storage  for  larger  M,  analytical  results  are  better  and  are 
the  recommended  approach  to  system  evaluation. 

The  investigation  here  has  assumed  that  only  one  sine  wave  is  present.  If  a 
second  sine  wave,  separated  in  frequency  by  more  than  the  width  of  an  NBF  in 
the  preprocessor,  is  present,  the  desired  tone  will  dominate  the  NBF  output  and 
the  current  results  will  apply.  Similarly,  if  the  preprocessor  consists  of  a  se¬ 
quence  of  FFTs,  separation  in  frequency  of  the  sine  waves  by  more  than  the  in¬ 
verse  of  the  time  duration  used  in  an  individual  FFT  will  not  cause  problems. 
However,  if  the  two  tones  lie  within  the  same  NBF  width  and  are  of  comparable 
strength,  phase-processing  can  suffer  severely  because  measurements  of  phase 
are  affected  by  both  sine  waves  and  will  vary  greatly  with  time.  In  this  situation, 
either  narrower  filters  should  be  used  or  amplitude  processing,  such  as  Proc¬ 
essor  III,  should  be  used.  This  situation  should  be  investigated  more  fully,  prob¬ 
ably  via  simulation. 


Figure  2A.  Receiver  Operating  Characteristic  for  PF  <  1 


Figure  2B.  Receiver  Operating  Characteristic  for  PF  <  -1 
Figure  2.  Receiver  Operating  Characteristic  for  Processor  II 
Known  Signal  Frequency 
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Figure  4A.  Receiver  Operating  Characteristic  for  Pp  <  1 
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Figure  4B.  Receiver  Operating  Characteristic  for  Pp£  • 1 
Figure  4.  Receiver  Operating  Characteristic  for  Processor  IV, 
Known  Signal  Frequency 
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Figure  5A.  Receiver  Operating  Characteristic  tor  Pp  <  1 


r 

Figure  5B.  Receiver  Operating  Characteristic  for  Pp  <  •  1 

Figure  5.  Receiver  Operating  Characteristic  for  Processor  V, 
Known  Signal  Frequency 
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C.  Receiver  Operatiikg  Characteristic 
for  c  (0,  2t r);  Band  (-*,*) 


3  6D.  Receiver  Operating  Characteristic 
t  (0,  2 7r) ;  Bands  (-*,*)  and  (0,2*) 

Receiver  Operating  Characteristic  for  Processor  VI, 
Known  Signal  Frequency 
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Figure  6E.  Receiver  Operating  Characteristic 
for  \p0  *  (0,  2tt);  Best  Tilted  Line 

Figure  6  (Cont’d).  Receiver  Operating  Characteristic  for  Processor  VI, 

Known  Signal  Frequency 
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Figure 


j  ha.  Receiver  Operating  Characteristic  tor  Pp£  1 


Figure  11B 


Receiver  Operating  Characteristic  for  Ppl  • 1 


non  11.  Receiver  Operatic  Characteristic  tor  Processor  II, 
Unknown  Signal  Frequency;  Track 
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Figure  12A.  Receiver  Operating  Characteristic  for  pFi  1 
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Figure  12B.  Receiver  Operating  Characteristic  for  Pp<  -1 


Figure  12.  Receiver  Operating  Characteristic  for  Processor  II, 
Unknown  Signal  Frequency;  Search 
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Figure  14A.  Receiver  Operating  Characteristic  for  Ppj£  1 
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Figure  14B.  Receiver  Operating  Characteristic  for  Pp  <  •  1 
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Unknown  Signal  Frequency;  Search 
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Figure  15C.  M  1 

Figure  15  (Cont'd).  Receiver  Operating  Characteristic  for 
Processor  VII,  Unknown  Signal  Frequency 
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Appendix  A 

DERIVATIONS  OF  PROCESSORS 


The  processors  will  be  discussed  in  the  same  order  in  the  appendixes  as 
they  were  in  the  text.  Also;  the  appendixes  are  arranged  such  that  particular 
processors  have  the  same  section  number  in  each;  e.g.,  Processor  V  for  known 
signal  frequency  is  treated  in  sections  A.  1.  5,  B.1.5,  and  C.  1.5. 

A.  1  KNOWN  SIGNAL  FREQUENCY 


A.  1.1  PROCESSOR  I:  KNOWN  SIGNAL 
PHASE;  AMPLITUDE  AND  PHASE  SAMPLES 

The  PDF  of  random  variables 

*  *[v  *mT,  &  >]T 

lor  signal  giesent  is  available  from  (3),  (4),  and  (5)  as 

1  I  6s, -T, CO/ ^4  Qk-P,  snri  i>aY 


p,U,y)  =  Tr  w 


2<rl 


(A—  1 ) 


(A- 2) 


The  PDF  for  signal  absent,  .  is  obtained  from  (A-2)  by  setting  TJ  =  0  . 

The  likelihood  ratio  (LR)  is  the  ratio  of  ]?,  to  p0  and  can,  using  (5),  be  put 
in  the  form 


lT  [  Km  J  2  r  J  • 


Therefore,  the  LR  test  is* 


uW*  JT, 


(A- 3) 


(A- 4) 


where  threshold  T  can  be  chosen  for  a  specified  false  alarm  probability.  Notice 
that  even  though  T0  is  assumed  known  in  (A-2),  it  is  not  used  in  the  LR  test 
(A-4). 


♦Satisfaction  of  the  upper  inequality  leads  to  the  decision  that  a  signal  is 
present;  satisfaction  of  the  lower  inequality  yields  the  decision  that  the  signal 
is  absent. 
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A.  1.2  PROCESSOR  II:  UNKNOWN  SIGNAL  PHASE; 
AMPLITUDE  AND  PHASE  SAMPLES 


When  the  signal  phase  is  unknown,  the  PDF  of  X  and  £  for  a  hypothesized 
signal  phase  *\>  is 

t>i  I  *■  -J 

If  the  a  priori  PDF  of  is  uniform,  there  follows 


Arr 


&7r  exr 


lx 

2«-'  J  0w'Llr 


(A- 6) 


The  PDF  for  signal  absent  is  obtained  by  setting  =  0  in  (A-6).  The  average 
LR  is  then  given  by 


(A-7) 


Therefore,  the  LR  test  is 


M  ;  Q 

2Xe  i  T 

K-l 


(A- 8) 


Again,  although  T0  is  assumed  known  in  (A-5),  it  is  not  used  in  the  LR  test. 

Instead  of  assuming  a  uniform  PDF  for  hypothesized  phase  ,  it  is  pos¬ 
sible  to  choose  'J'  as  the  ML  estimate  in  (A-5).  By  rearranging  (A-5)  intheform 

T,  U ,  alt) «*[-  +  |rTi.[e j  >  ,A-9) 

it  is  readily  apparent  that  p,  is  maximized  by  the  ML  estimate 

$>  =  e1*}  .  <A-10) 


(A-10) 


The  maximum  value  of  p,  is  then  given  by 

,  i  r  v 


F 

^  )c-t  J  ■ 


(A —11) 
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The  generalized  LR  is  the  ratio  of  (A- 11)  to  pe  : 

e*p["  2?^  1]  • 


(A — 1  2) 


The  generalized  LR  test  is  therefore  identical  to  (A-8).  Thus  two  different 
methods  of  treating  signal  phase  yield  the  same  processor  (A-8). 


-.1.3  PROCESSOR  III:  INDEPENDENT  SIGNAL 
PHASES;  AMPLITUDE  AND  PHASE  SAMPLES 


When  all  signal  phases  are  independent  from  sample  to  sample,  the  PDF  of 
x  and  Q  is  given  by  a  modified  form  of  (A-2),  namely, 


i  ^  -1 


(A- 13) 


where  is  the  set  of  M  random  signal  phases.  For  uniform  PDFs 


of  each  v^K  ,  there  follows 


p  (x,  a)  =  (i^f  t.U,  i  li) 

,  -ft  M 

- _ ! _ _  eyp  &**  +  *?•  TTt  (^~t) 


(A —14) 


The  average  LR  follows  readily  from  (A- 14)  as 
,  M  , _ 


«P(- 


(A -15) 


and  the  LR  test  becomes 


i  T 


(A- 16) 


This  test  makes  no  use  of  phase  samples  ,  but  does  depend  on  knowledge  of 
^  ;  however,  for  small  SNR  (small  T/cr  ),  (A -16)  becomes 

JL-Rk  >  T  (A-17) 


For  small  SNR,  this  approximate  LR  lest  does  not  require  knowledge  of  T0  . 
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A.  1.4  PROCESSOR  IV:  KNOWN  SIGNAL  PHASE; 

PHASE  SAMPLES 

The  PDF  of  amplitude  and  phase  samples  and  is  available  from  one 
term  of  (A-l)  by  transforming  to  polar  coordinates: 


?.(Mk>  «p[-  ],VA  M 

The  PDF  of  each  9^  is  obtained  by  integrating  over  "Rs  : 


<TT.  (A-l  8) 


?,  f©K>  =  j^[«p(-  ^r)  +  J-  no(eK-lt.).’xp(-S.s',„<(eK-'ll^ 


“I  , 

J  4*  exp(-V/*)  ;  K\  <  ^ 


(A— 19) 


Now,  letting  Ch-Co^&K-^  for  notational  simplicity,  we  expand 

?,  (6.)  -  ^  [l  +  Vf  f  G  +  i)  +  O  (-^r)]  ,  kP-TT,  (A-20) 

where  0(  )  denotes  terms  the  order  of  (  )  . 

On  the  other  hand,  consider  the  approximation  to  p,  (O  of  the  form 

f,W  -  (A-21) 

This  approximation  has  unit  area  and  the  expansion 

M»k)  -  |-G  +  ?i(C-i)+0($.)],  W<ir.  (A- 22) 

Now,  (A-20)  and  (A-22)  are  identical  through  order  T^/r,  and  are  almost  equal 
through  order  T^/V  ,  i he  factor  1  in  (A-20)  being  replaced  by  ttA-  in  (A-22). 
Therefore,  we  use  the  approximation  (A-21)  from  this  point  on.  Madiematically, 
we  write 

f,  (O  *  **y  fit  A(vf  ^  orcJer  ^  (A-23) 

while  remembering  that  (A-23)  is  also  a  good  approximation  to  order  T0  / ^  . 
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The  joint  PDF  of  [aj  is  then  given  by 


?,(e)=^rexp^s|c,]/^(1|'i)j4  o^«r  £,  I6kHx 


(A -24) 


The  PDF  of  jj)  for  signal  absent  follows  from  (A-24)  by  setting  T0  =  0  .  The  LR, 
to  order  V0/v~  ,  is  therefore  given  by  (A-24),  with  the  factor  absent.  Thus, 

the  LR  test  is 


JElC  =  R«[e  'Vo  iEe'6*]  <  T,  4  order 

i  *  L  k* i  J 


(A-25) 


and  is  independent  of  “P#  and  <r .  Therefore,  threshold  T  can  be  selected  once 
and  for  all  to  realize  a  specified  false  alarm  probability,  ~PF  .  The  value  of  T 
depends  on  only  M  and  TF  ,  and  Processor  IV  is  a  C FAR  receiver;  i.e.,  test 
(A-25)  is  uniformly  most  powerful  wi  h  respect  to  and  or  ,  to  order  "P0/<r~  . 


A. 1.5  PROCESSOR  V;  UNKNOWN  SIGNAL 
PHASE;  PHASE  SAMPLES 


The  PDF  of  0  for  unknown  signal  phase  is  given  bv  (A-24),  where  in 
Q<  is  replaced  by  hypothesized  angle  ;  i.e. , 


f.(sH')  =  ^  ^r[ff  /tT(V?' £-), 

•to  Order  >  |9K|<"7r- 


(A-2f>) 


We  now  express 


ReL'^ieS  . 

jj-i  (.  K1 1  9 


(A-27) 


Then,  for  an  a  priori  PDF  of  phase  \|)  that  is  uniform  over  2tt,  the  PDF  p,  (§) 
is  given  by 


T»(S) 


2tr 


2-n- 


(A-28) 
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Since  the  PDF  of  6  for  signal  absent  is  obtained  from  (A-28)  by  setting  ^>  =  0  , 
the  average  LR  is  given  by  (A-28)  without  the  factor  (2-rr)-M  •  Therefore  the  LR 
test  is 


S  ^  ^  T(  4d  ortJer  -jp- 

fc*  i 


(A-29) 


If,  instead  of  assuming  a  uniform  PDF  for  >)j  ,  we  choose  ijj  as  the  ML  es¬ 
timate  in  (A-26),  we  find,  using  (A-27),  that  the  ML  estimate  is 


a  y  *  ie*) 

+  ’  a"3l^e  1 

and  the  maximum  value  of  "p,  is 

T.ffil  V)  *  ^  e*r  i  e,e"  |]  /t;  (VF£) , 

-(«>  o  r  cW  >  1 0*  I  ^  TJ“, 


(A-30) 


(A-31) 


The  generalized  LR  is  the  ratio  of  (A-31)  to  Qbr)  and,  again,  the  generalized 
LR  test  is  obviously  (A-29).  Thus,  to  order  P„/cr  ,  two  different  methods  of 
treating  signal  phase  yield  the  same  processor.  Test  (A-29)  can  be  designed  for 
specified  Yf  without  knowledge  of  or  <r. 


A.  1. 6  PROCESSOR  VI:  UNKNOWN  SIGNAL 
PHASE;  FITTED  PHASE  SAMPLES 

Instead  of  attempting  to  derive  average  LR  or  ML  tests,  we  consider  here 
a  heuristic  approach  to  detection.  Specifically,  th 2  M  successive  [0^]  are 
fitted  by  the  best  straight  line  such  that  the  average  squared-error  is  minimized. 
Then  the  resultant  minimum  error,  or  scatter,  is  used  as  a  decision  variable. 
For  small  SNR,  the  phase  samples  would  be  widely  scattered  over  a  2-ir  inter¬ 
val,  whereas  they  would  tend  to  cluster  around  the  true  signal  phase  for  large 
SNR,  yielding  little  scatter.  Thus  if  the  scatter  is  less  than  a  threshold,  we  de¬ 
clare  the  signal  present. 

For  known  frequency,  the  best  straight  line  has  zero  slope,  as  can  be  seen 
in  (5).  Therefore  the  average  squared-error  for  hypothesized  phase  'p  is 


This  is  minimized  by  the  choice 


(A-32) 


(A-33) 
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which  is  the  average  phase.  The  minimum  error  is  therefore 


(A-34) 


A 

When  E<  T,  the  signal  is  declared  present.  Test  statistic  (A-34)  is  independent 
of  and  <r  and  yields  a  CFAR  processor;  the  samples  need  not  be  equally 
spaced  in  time. 


There  is  a  problem  in  test  (A-34)  as  to  deciding  the  2v  interval  or  band  in 
which  to  choose  each  ©K  .  The  choice  definitely  affects  detectability,  the  exact 
degree  depend'ng  on  the  true  (unknown)  signal  phase  .  This  problem  is  dis¬ 
cussed  fully  in  the  main  text. 


A.  1.  7  PROCESSOR  VII:  UNKNOWN  SIGNAL 
PHASE;  PHASE-DIFFERENCE  SAMPLES 

No  derivation  is  necessary  for  this  case.  The  phase-difference  samples 

are  formed  for  2sKarL  For  known  signal  frequency,  the  phase  differences 
tend  to  cluster  around  zero  for  high  SNR,  even  though  signal  phase  < pp  is  unknown. 
Thus,  a  random  walk  similar  to  (A- 25)  is  formed  and  the  real  part  is  compared 
with  a  threshold.  The  test  is 

exj>[i  (0K-  £  T.  'A -35) 

The  samples  need  not  be  equi-spaced  in  time. 

A.  2  UNKNOWN  SIGNAL  FREQUENCY 

It  is  unrealistic  to  assume  that  signal  phase  would  be  known  at  some  time 
instant  for  unknown  signal  frequency.  Accordingly,  there  are  no  analogs  of 
Processors  I  and  IV  in  this  case. 

The  received  waveform  here  is  a  slight  modification  of  (3): 

r  It;  -  "Re  [  exp(i  2rr(f-0-'-  \Y0  (i  ^0)  +  HI  (£)]]  ( A  -  3  6 ) 

for  Doppler  shift  .  The  complex  envelope  of  the  received  process  is 

exp  (i2Trfdt)[?0  exp(li)  +  13  (,*)]  (A-37) 

A  sample  of  the  complex  envelope  at  time  -tM  is  denoted  by 
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\e  "  “  ^H4'^  =  e*F (i^irf,,**)^  ftXfOk)  +  J3  fe)J 

(A-38) 

-v-  n„  u^i^fo. 

The  statistics  of  the  real  component  samples,  nxftM)  a^d  O^tty,)  ,  of  the  noise 
are  available  from  the  following  (reference  3): 

\{k)  +i^(t)  =  «'fc)'  txp(i2irfj-t)  =  0, 


)  ft)  +  i  ^  f£)|  = 

(Vlx  (fc)  +  i  rtj  rt)]7  =  £  ft)  ex^fi^-n-fj-t)  *  0. 

Therefore,  _  _ 

r>;ie  -  *<ji9  =  o-, 

_ _ _  (A-40) 

=  0. 

These  noise  statistics  are  identical  to  those  for  known  signal  frequency. 


Inter  -  2<r* 


(A-39) 


A.  2.1  PROCESSOR  I:  NO  ANALOG 


A.  2.  2  PROCESSOR  II:  UNKNOWN  SIGNAL 
PHASE;  AMPLITUDE  AND  PHASE  SAMPLES 


For  signal  present  with  hypothesized  phase  i{j  and  frequency  shift  -f  ,  the 
PDF  of  samples  \  and  ^  in  (A-38)  is,  using  (A-38)  and  (A-40), 


P'(X*  yJt,-f)=  2  ~~~T  T-p  ^  ^VT° ** 


(A-41) 


The  PDF  of  x  and  ^  in  (A-l)  is  therefore 

M*,#,  f)  '  6t[- 


M 


f A —421 
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(A-43) 


If  the  a  priori  PDF  of  t|t  is  uniform,  there  follows 

If,  at  this  point  we  make  an  assumption  about  a  PDF  for  -f  ,  we  must  make  some 
approximation  to  (A-43)  so  that  the  evaluation  will  be  amenable.  The  one  weadopi 
is  identical  to  that  made  earlier,  namely,  small  SNR.  Then,  the  relevant  inte¬ 
gral  hat  must  be  evaluated  is 

cK  «  *p[|  M  ^  ^  4  order  ^0  >  (A_44) 

where  & ,  +»)  is  the  range  of  anticipated  Doppler  shifts,  and  the  PDF  is  assumed 
to  be  uniform.  Now,  fora  large  range  f, -f,  ,  the  only  terms  in  (A-44)  that 
yield  significant  contributions  are  K=JJ  : 


4  order  * 


(A-45) 


Alternately,  (A-44)  again  reduces  to  (A-45)  for  equi-spaced  samples, =4  , 
and  a  uniform  PDF  of  f  over  a  '/a  Hz  range  (which  is  approximately  the  fre¬ 
quency  separation  between  adjacent  filters  as  discussed  following  (2)).  Thus  no 
use  would  be  made  of  [0^  under  either  assumption  for  the  PDF  of  f  .  This  is 
a  case  of  designing  for  the  worst  possible  situation  and,  thereby,  drastically 
degrading  performance  (see  reference  3,  page  180). 

A  better  t^proach  is  to  return  to  (A-43)  and  choose  -f  as  the  ML  estimate; 
i.e.,  choose  f  such  that 


i6s 


>  Ke  e 


allf*  MO,  (A"46) 


where  is  the  allowed  range  of  Doppler  shifts.  Then  the  comparison  of 

I  %<*,*)  with  a  threshold  is  equivalent  to  the  test 

>  T.  (A-47) 

K=  i 

Another  approach  is  to  use  ML  estimates  for  both  vjj  and  f  .  From  (A-42) 
we  see  that  the  ML  estimate  of  vjj  is 


;  f  10x 
t  =  e. 


(A-48) 
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Substitution  of  (A-48)  in  (A-42)  and  the  subsequent  choice  of  f  for  a  maximum 
of  7,  again  yields  (A-46);  and  the  resultant  generalized  LR  test  is  age- in  (A-47). 
This  test  is  independent  of  ;  also,  no  assumption  about  small  SNR  is  required 
to  deduce  (A-47)  as  an  appropriate  test. 


A.  2.  3  PROCESSOR  III:  INDEPENDENT  SIGNAL 
PHASES;  AMPLITUDE  AND  PHASE  SAMPLES 


For  independent  signal  phases,  the  PDF  is  a  modified  form  of  (A-41): 
fi  (\,%  "0  ~  2^-<rx  e'Xf‘jj  20^1^*'  Cos(2nA{K+  ^ 


28-* 


[yK-  ^  sin 


(A— 49) 


Therefore, 


(A-50) 


(A—  51 ) 


U 

Now,  if  the  phases  are  uniformly  distributed  over  2-rr , 

f.  u  *  13  ■  -r[-  M  • 

Since  this  PDF  is  independent  of  -f  ,  a  decent  processor  can  not  possibly  result 
and  we  discard  to  is  approach. 

If  we  instead  take  the  ML  estimates  for  ^  ,  we  find  from  (A-50)  that 

k  (A- 5 2) 


Then 


fiUsl*,3=  exr 


M 


(A-53) 

[  2s-  tc)  2r  <r  k* ,  J. 

Again,  the  independence  of  f  yields  an  undesirable  processor.  Notice  from 
(A-53)  that  the  generalized  LR  test  would  take  the  form 

i  ^  X  (A-54) 

t 
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A.  2.4  PROCESSOR  IV:  NO  ANALOG 

A.  2.  5  PROCESSOR  V:  UNKNOWN  SIGNAL 
PHASE;  PHASE  SAMPLES 

Our  starting  point  is  (A-41).  By  transforming  to  polar  coordinates  according 
to  (5),  we  obtain 


Then 


?,kM ■  4  onW  *-'  ,A_56) 

by  means  of  an  approach  that  is  analogous  to  that  in  (A-18)  through  <A-a>. 

Th6ref0re  fr _ _  "P  'e«  -‘^11  x 

7.(il^)  =  K--ie  e  ^  (A-57) 

If  we  assume  that  the  PDF  for  4»  is  uniform,  we  obtain 

r.(al-f)  =^L(i?  £-1  i^rOe«-'2-M) *  <*-58> 

x;M  (i/?  ,  -t*  •***■  5-  • 

Accordingly,  we  choose  the  ML  estimate  -f  m  (A-58): 

\±?«M-^\>-\  t^*K i9«-;^L)l.  •»  *'U\  <A-59) 

and  the  generalized  LR  test  then  takes  the  form 


(A-58) 


8W  -  i  2-tt'F^:k)|  <  T)  4&  otrdefc'  T' 
f,<f<K  k=i  ' 


(A-60) 


The  alternate  approach  of  using  ML  estimates  for  both  .V  and  f  in  <A-57> 
is  also  possible:  we  find  that 

I  ,  arg  [  Jtexr(ieK-i  2rfO]  A'*'-?-'  (A'61) 
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Substitution  of  (A-61)  in  (A-57)  and  the  subsequent  choice  of  -f  again  leads  to 
(A-59)  and  the  generalized  LR  test  (A-60).  This  test  is  independent  of  T0  ,  but 
is  accurate  only  to  order  1^/tr  ■ 

A.  2.  6  PROCESSOR  VI:  UNKNOWN  SIGNAL 
PHASE;  FITTED  PHASE  SAMPLES 


The  philosophy  here  has  been  explained  earlier  in  this  appendix.  However, 
the  average  squared-error  for  unknown  frequency  is  generalized  from  (A-32)  to 


E  * 


(A-62) 


W  K-  i 

where  we  assume  that  the  phase  samples  are  equi-spaced  in  time.  The  partial 
derivatives  of  E  with  respect  to  vjy  and  jS  are  both  set  equal  tc  zero  and  solved 


to  obtain  the  estimates 


K=  I 
6 


K=,  •<  m(m-o  rr, 

Substitution  of  the  estimates  (A-63)  in  (A-62)  yield  the  minimum  scatter 


(A-63) 


A 


E  = 


(A-64) 


When  £  <T,  the  signal  is  declared  present.  Test  (A-64)  is  independent  of  Ta 
and  is  a  CFAR  processor. 


A.  2.  7  PROCESSOR  VII:  UNKNOWN  SIGNAL 
PHASE;  PHASE-DIFFERENCE  SAMPLES 


For  unknown  signal  frequency  and  equi-spaced  samples,  the  phase  differences 
0K-&K.1  cluster  around  an  unknown  angle  for  large  SNR.  T  le  value  of  the  unknown 
angle  depends  on  the  signal  frequency  and  the  time  betwe<  a  samples.  A  random 
walk  and  magnitude  similar  to  (A-29)  is  formed  since  only  the  length,  and  not 
the  direction,  indicates  signal  presence.  Thus  the  test  is 


$  T. 


Again,  the  test  is  independent  of  and  yields  a  CFAR  processor. 


(A-65) 
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Appendix  B 

GEOMETRICAL  INTERPRETATIONS  OF  PROCESSORS 


No  assumptions  about  Gaussian  noise  or  independent  samples  are  made  in 
this  appendix.  Rather,  geometrical  interpretations  of  "good"  processors  are  de¬ 
veloped  and  the  errors  associated  with  each  are  minimized.  The  resultant  tests 
are  identical  to  those  derived  in  appendix  A.  Although  signal  amplitude  is 
assumed  known,  it  is  never  needed  in  the  tests. 


B.  1  KNOWN  SIGNAL  FREQUENCY 

B.  1.  1  PROCESSOR  I:  KNOWN  SIGNAL  PHASE; 
AMPLITUDE  AND  PHASE  SAMPLES 


(  Uy 

If  a  signal  is  present,  the  complex  samples  1?he  should  cluster  about  the 
point  ?0  e1*0  .  A  measure  of  the  scatter  is  afforded  by  the  sum  of  the  squared  dis¬ 
tances  between  these  points: 

e,» 

---  It  ^M-P„’-2R  T?e[ei<4  ibA”']  .  <B  ^ 

If,  on  the  other  hand,  the  signal  is  absent,  the  complex  samples  should  clus¬ 
ter  about  the  point  0.  The  scatter  is  then 

_  4r-  U  »  r-,1  (B-2) 


E„  .  X  |Ke- 


-0|  = 

k=  i 


Now,  if  E,  <  E0  by  a  sufficient  amount,  we  would  be  quite  sure  that  a  signal  is 
present.  If,  on  the  other  hand,  the  converse  is  true,  we  would  declare  that  the 
signal  is  absent.  Therefore,  the  test  we  adopt  is* 

£>  >  E0  ~  A.  (B-3) 

Substitution  of  (B-l)  and  (B-2)  in  (B-3)  yields 

1?e$Le  Jt-KX*"}  5  T,  (B-4, 

where  T  is  a  threshold  that  is  adjusted  for  a  prescribed  PF  .  Test  (B-4)  is 
identical  to  (A-4). 

*Here,  as  in  the  remainder  of  this  appendix,  satisfaction  of  the  upper 
inequality  leads  to  a  statement  of  signal  presence,  whereas  satisfaction  of  the 
lower  inequality  leads  to  a  signal-absent  decision. 
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B.  1.  2  FROCESSOR  II:  UNKNOWN  SIGNAL  PHASE; 
AMPLITUDE  AND  PHASE  SAMPLES 

For  a  hypothesized  signal  phase,  ^  ,  the  scatter 

_  iu  ie-  «  "VI2 


E,  =  ^  V  -?.e 


(B-5) 


We  can  select  so  that  tie  scatter  is  minimized;  i.  e. ,  we  try  to  find  that  com¬ 
plex  point  l^exiM)  about  which  the  samples  cluster  best  if  a  signal  is  present. 

The  value  of  that  minimizes  (B-5)  is 


(B-6) 


(Note  that  this  is  identical  to  the  ML  estimate  A-10. )  Then 

e, mt£  -  zr.  I  itiuiei 


(B— 7) 


For  signal  absent,  E„  is  again  given  by  (B-2).  Substitution  in  (B-3)  yields 


ih,e:e* 


>  -p 
<  '  > 


(B-3) 


which  is  identical  to  (A-8). 


B.  1.  3  PROCESSOR  IH:  INDEPENDENT  SIGNAL 
PHASES;  AMPLITUDE  AND  PHASE  SAMPLES 

For  hypothesized  6  gnal  phase  in  sample  k  ,  the  appropriate  measure 
of  .scatter  for  signal  present  is 

it  2 

E,  =  -£  I?**  -Fo« 

_  . ,  .  M  (B-9) 

=  XK-vM?:-2P„'R.^Ke'  *'“] 

Since  we  must  be  allowed  to  minimize  E,  by  choice  of  \JiK  ,  we  choose 

A 

(B-10) 

yielding  M  M 

E,  =■■  +  M  ?l- 2P„  5.P* . 
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If-  t 


*  f 


I  l 

it.  t 


Scatter  is  given  by  (B-2).  Test  (B-3)  then  yields 


K-i 


K 


>  -T- 
<■  »  * 


(B-12) 


Test  (B-12)  is  not  identical  to  LR  test  (A- 16)  or  to  the  approximate  LR  test  (A- 17) 
that  was  derived  lor  small  SNR,  However,  if  the  SNR  is  large  i)  ,  (A-16) 

tends  to  (B-12;  and,  therefore,  (B-12)  is  nearly  optimum  for  high  SNR.  Test 
(B-12)  is  not,  however,  pursued  further  in  this  report. 

B.  1.4  PROCESSOR  IV:  KNOWN  SIGNAL  PHASE; 

PHASE  SAMPLES 

Since  amplitude  samplt.1'  are  not  used,  it  is  expected  that  complex  sam¬ 
ples  would  cluster  around  the  point  for  large  SNR.  Therefore,  a 

meaningful  measure  of  scatter  is  given  by 

M 

Ei  = 


K  -  I 


e  -  e 


Scatter  En  is  now  given  by 


E.  ~  I  £  -  0 


■- 1 


'0-K 


=  M. 


Thus,  evaluation  of  (B-3)  yields 


-'A 


e )  e 


^  T> 


which  is  identical  to  (A- 25). 


B.  1.  5  PROCESSOR  V:  UNKNOWN  SIGNAL 
PHASE;  PHASE  SAMPLES 

For  hypothesized  signal  phase  \Ji,  the  scatter  for  signal  present  is 


E.  - 


M 


K=  i 


.  e  -  e. 


"V 


=  2M 


K-  i 


This  is  minimized  by  the  choice 


(p  =  Q^rq 


io, 


K-  I 


(B-13) 


(B-14) 


(B-15) 


(B— 16) 


(B-17) 
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which  is  identical  to  ML  estimate  (A-30),  and  leads  to 

E,'2M-2|±eifct 

Scatter  £  is  given  by  (B-14)  and  substitution  in  (E-3)  yields 


n 


\B, 


\ 


>  T, 


which  is  identical  to  (A-29). 


(B-18) 


(B— 19) 


B.  2  UNKNOWN  SIGNAL  FREQUENCY 

It  is  unrealistic  to  assume  that  signal  phase  would  he  known  at  some  time 
instant  for  unknown  signal  frequency.  Accordingly,  there  are  no  analogs  of 
Processors  I  and  IV  in  this  case. 

B.  2.  1  PROCESSOR  I:  NO  ANALOG 


B.  2.2  PROCESSOR  II:  UNKNOWN  SIGNAL  PHASE; 

AMPLITUDE  AND  PHASE  SAMPLES 

For  high  SNR,  the  complex  samples e'^]  in  (5)  should  cluster  around  a 
uniformly  rotating  complex  vector,  P0  y  \  Trrft)  ,  where  initial  phase  ^  and 
rotation  rate  -f  are  unknown.  The  measure  of  scatter  adopted  is  therefore 


E,  =  jS-  Ne’6"  -  k  exp  (if+ 

-  +  M?’  -  2k  p4e  4  21V'e’  e 

k=i  0  l  K-i  K 


t  TiT-f-tKj 


(B-20) 


This  is  minimized  by  the  choice  of  hypothesized  signal  phase  as 


h  aT 


XR, 


;eK  -  i  >rrfO 


(B-21) 


which  is  identical  to  (A-48).  The  resultant  value  of  E.  , 


M 


E,=  ^K  +  M?J-2V0 


K-i 


M. 


\9*  - 1  TrrTF^ 


(B-22) 


is  further  minimized  by  choosing  f  such  that  the  last  term  in  (B-22)  is  maxi- 

A 

mized;  i.e.  ,  choose  f  such  that 
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,6,  -'.irf^l  |«  ift.  -IW-i,  ,.  .  ,.  . 

£  j  -  vf  6  -for  all  Hr  f  (■pj'R.)-  (B-23) 


This  is  identical  to  (A-46).  Then 

E,  =  2.^  +  MP0-2f?  «“*  .SV  & 

K - 1  "  k-i 


(B-24) 


The  value  of  scatter  for  signal  absent  is  still  given  by  (B-2).  Therefore, 
(B-3)  yields 

1 0K  - 1  2.Trf 

rw  2-\z  e  >  T  (B- 

k»i  K  ’ 

which  is  identical  to  (A— 47). 

B.  2.  3  PROCESSOR  III:  INDEPENDENT  SIGNAL 
PHASES;  AMPLITUDE  AND  PHASE  SAMPLES 


(B-25) 


The  approprate  measure  of  spread  for  independent  phases  is 
i—  ^rlo  -X  ^  Xi7 


r  ^T|o  X  ~  X)7 

E,  =  S|K„e  -?„e  |  _ 


(B-2  G) 


Since  this  is  the  same  measure  that  was  used  for  known  frequency  in  (B-9),  test 
(B-12)  is  not  an  acceptable  test  for  unknown  frequency;  it  is,  however,  identical 
to  (A- 54), 

B.  2.  4  PROCESSOR  IV:  NO  ANALOG 

B.  2.  5  PROCESSOR  V:  UNKNOWN  SIGNAL 
PHASE;  PHASE  SAMPLES 


The  appropriate  scatter  here  is  a  modification  of  (B-20): 
*1  (0* 


C  Jj  -i'iirfi 

=  2M-2Keey-^  e 


(B-27) 
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M.  scatter  is  minimized  by  the  chotce  of  hypothesized  signal  phase 

$  =  (B'2S) 
vhich  is  identical  to  (A-61).  Then 

E,  =  2M-2^e*p(jeA'*TfO| 

can  be  further  minimized  by  choosing  {  such  that 

|  1_  ei6“  >  Hie1®'  4r  .11  4  £  (VA 

which  is  identical  to  (A-59).  Then 

E,  =  2  M  -  2  "ft  1  %  «t0  *•  - 1 2irft:)l  • 

The  vaiue  of  scatter  for  signal  absent  is  still  given  by  (B-14).  Therefore 

3)  yie,ds  \  m  .  \ 

2  0K?(l'6a- T 

f- 1 


(B-29) 


(B-30) 


(B-31) 


vn«v 


(B-32) 


ch  is  identical  to  (A-60). 
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Appendix  C 

DERIVATIONS  OF  PROCESSOR  PERFORMANCE 

The  additive  noise  is  assumed  to  be  Gaussian,  and  the  samples  statistically 
independent  as  noted  in  the  discussion  following  (2).  Derivations  of  the  detection 
probability,  ,  will  be  made  or  approximated  where  possible,  and  numerical 
methods  that  would  lead  to  exact  performance  calculations  in  some  cases  will  be 
outlined.  Such  situations  were  not  pursued  to  completion  in  all  cases,  however. 

Some  of  the  derivations  to  follow  are  condensed  in  the  interest  of  brevity. 
However,  the  essential  steps  are  presented  in  such  a  way  that  the  reader  can 
follow  them  and  insert  the  missing  calculations. 

C.  1  KNOWN' SIGNAL  FREQUENCY 

C.  1.  I  PROCESSOR  1:  KNOWN  SIGNAL  PHASE; 

AMPLITUDE  AND  PHASE  SAMPLES 


From  table  1  and  equation  (5),  we  express  the  decision  variable  as 

M  M 


ri  —  i  k  r  i 


(C-l) 


( C—  2 ) 


Random  variables  (RV)  £  and  are  Gaussian,  with  the  PDF  given  by  (A-2) 
for  signal  present.  Then  RV  X  is  Gaussian,  with 

£  (JH)  =  cos  v[>o  M?„  cos<|>p+ StVnJi  MV0  s  im  ^  —  Mft>, 

Var^  =  c*s>0  N\  Varfv)  d  M  Var(x)  -  M 

where  we  used  the  independence  of  all  the  RVs.  Therefore 

PD  -  Pr.tD>T)  = 


(C-3) 


where 


x  r  y 

$(x)  s  J  i,W  t<ph'/’)u3  W.  (C_4> 

—  °o 
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T  .  .  % 


If  we  define  a  normalized  threshold  X  ~  and  se*-  H  follows  from 

(C-3)  that  1 


T.-K.t/S’wi-v),  PF  *!(->>); 

false  alarm  probability  PP  was  obtained  by  setting  w\~o  in  PD. 

C.  1.  2  PROCESSOR  II:  UNKNOWN  SIGNAL  PHASE; 
AMPLITUDE  AND  PHASE  SAMPLES 

From  table  1  and  (5), 

i-  [(£4  •  &•)']*. 


(C-5) 


(C-0) 


The  PDF  of  X  and  £  is  given  by  (A-2),  and  since  U.  and  V  are  Gaussian  RVs, 
we  have 

ProW  U  >  ^  =  ^  ^U/ ^ 


=  tadv  e-f  -  ^p[(lA-M^C05^o)2+ 

{[ 7+?>t 

0° 


1 


=  f  drr  [de^TMrJl«xf»[~a?r^+d^--2MP0rcaS(6^ 


,  r  I  r’+h’lflx  CEhP} 

^  FT"  exrlr  _2FF^JI°  CF7 


( C— 7) 


T 


-  q(^>^  *  G‘(vffm> 


independent  of  ,  where  (reference  3,  appendix  F) 

Q(qi  b)  =  j^dx  x  exp(-  x  J  "  )  T0  («k).  (C-H) 

The  re  fore 

?D  =  Q  (i/m  *»,  M  PF  =  Q  (0,  \)  -  exP(-  >>%).  <c-9» 


■p  ^ '  jt^wvTjTji  «  r^fYK'ir:w7ft"r*T jjjj 
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C.  1.  3  PROCESSOR  III:  INDEPENDENT  SIGNAL 
PHASES;  AMPLITUDE  AND  PHASE  SAMPLES 

From  table  1  and  (5), 

M 

(C-10) 

K  =  ! 

The  characteristic  function  (CF)  of  RV^(  is 

f(?)=  =  TT  E[e*p[i?(xK2+y*)]j 


upon  transforming  to  polar  coordinates  and  using  reference  (i,  li.ii.31  4.  The 
probabilities  of  detection  and  false  alarm  are  then  (relerenee  3,  pp.  217-219) 


"  Qm  {Jt?  X4  ;  \)  ,  '  /*  (  °)  > 

where 

''»v(  "V")  »„  ,H 

-  Q(#,t)+  «P(  «';>'■)  u«w, 


(CM2) 


(C- 13) 


and,  in  particular, 


%  =  ^  =  exi'(  M  o  ^  ; 


Extensive  numerical  results  are  provided  in  releroin 


C.  1.4  PROCESSOR  IV:  KNOWN  SIGNAL  I’llASE; 
PHASE  SAMPLES 


(C-14, 


It  is  convenient  at  this  point  to  make  a  change  of  variables;  using  (:>), 

(x„  ■+  i  y«)  exp(-  Ife)  =  R^xp(ie.-d.)  -  ?.  +fl(t> ,<4  E  TJ  4-<r(a„+i  k)(C'lr” 
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That  is, 

aH  +  i  K  =  r"  Jrt  fe)  ■  (C-16) 

Then  (reference  2  or  4) 

E  (a*-)  +  '  E  (l*')  =  Jr  exp  (-1  >li)  E  (  a  Ik))  =  0, 

E<|M>Wr)=  EW  +E(W) ->E(|»(tjr)  --  2,  <c-17> 

E(.^+ib,7)  =F(ce)E(0+i2E(\K)  =  ^e'2)lE(ttfe))-  0. 


Therefore 


Ek)=  E(W)  =  o, 

eKW)  =  o, 
e(«;)=  e(w)=  '• 


(C-1H) 


Since  (C-1G)  is  a  linear  transformation  on  a  complex  Gaussian  RV,  (C  18)  indi 
cates  that  QK  and  kK  are  independent  zero-mean  unit-variance  Gaussian  ItVs. 


Then  from  table  1  and  (C-15),  the  decision  variable  becomes 


1  4P 

hi  fox 

If 

_ 

[W+kfl*  ' 

where 

Exact  Detection  Probability 
The  first  order  PDF  of 
S  h 


VV>  H 


T 
1  0 

<r 


(C- 19) 


(C-20) 


F  bl]  ^ 

can  be  evaluated  in  closed  form  as  follows.  Let  t=  ;  then,  for  lrl< 

Pro  \>  C3  <  r)  =  Pro^  (  <  *)  =  ^rol> 


(C-21) 


-iu) 


(C-22) 


_ «  -  *0  '  -M  -*0  -*> 


(i  1 


si-.*  ■  *7  i®5<wa^is*^Ews?!  y»*^'i^3^g»E^^y^^^<?g^ywfe9^5>'**igigg;^r??,i?;?g^wg^iBggg¥gg>»i^gsg|^ 
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using  (C-21),  (C-18),  ;tnd  (C-4).  Therefore,  the  PDF  of  >S,  for  |r-|<^  I,  is 

fc(0  -  ^rfrol(s<  v)  =  _j<H>  4>0>)  5 '  [j==F  ~  w)  -t  jyj^} 

.(^rj>w(js!Mw 


(C-23) 


==:  txpf-kyiTi)  4-/2tt  wirexpf-(i-r*)wy^  3>(*ir) 

-r  1  v  '  v  _  > 


Vi-r1 


upon  completing  the  square  in  the  exponent. 

Now  the  CF  of  S  ,  which  is  given  by 

fs  ff)  =  /^expOf^psfr),  (C- 

can  be  quickly  evaluated  by  a  fast  Fourier  transform  (FFT).  (This  CF  is  also 
available  directly  from  (C-21)  and  (C-18)  as 


(C-24) 


•fjff )  *  f°dr  r  exp  (-  r  -  j;^1-)  I.  (vnr  +•  i ?) , 


(C-25) 


but  that  method  is  much  more  time-consuming.)  Then,  the  CF  of  RVJJ  in  (C-19) 
is  -fsM(r).  Finally,  the  cumulative  distribution  of  RV^  is  available  (from  ref¬ 
erence  7,  (7))  as 

PD  =  TVok(i>-r)  -  i  +  ^(^Tm[-f"(f)exp(-ifT^  <c-26> 

which  c;m  also  be  evaluated  by  an  FFT  (reference  8).  This  exact  pioceduxe  has 
not  been  pursued  in  this  report. 

Exact  False  Alarm  Probability 


The  CF  of  s  for  signal  absent  is  given  by  (C-25)  with  tvi  =  0  ■ 

f,(?)=X„(i?)  =  T0(?). 

Therefore  the  CF  of  B.vJ  in  (C-19)  is 

-gfr)  =  o^(r). 


(C-27) 


(C-28) 
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and  thus  (reference  7,  (7)),  v, 

TF  .  ?,-0^>T)  -  i  -  MC  ft®  ^Irr) 


where  normalized  threshold 

\  S 


_M 

0 


(C-29) 


T. 


(C-30) 


Equation  (C-29)  is  tin  exact  relation  for  ?r  ;  it  is  seen  to  depend  on  only  X 
and  M  Therefore,  Processor  IV  is  a  CFAR  receiver;  i.  e. ,  X  can  be  chose 
to  realize  a  prescribed  PF  ,  requiring  only  knowledge  of  M  .  The  reason 
the  threshold  scaling  in  (C-29)  and  (C-30)  is  that 

Thus,  (C-29)  yields  (see  referenced,  3.896  4;  integrate  both  sides  with  respect  to  b) 

F>~  <£(-*)  «  (c~32> 

Therefore  it  is  anticipated  that  the  choice  of  X  in  the  exact  relation  (C-29)  for 
specified  PF  becomes  relatively  independent  of  M,  for  large  M  .  table  C 
gives  values  of  X,  as  determined  from  (C-29),  for  =  25"  . 

Table  C-l.  Thresholds  Required 
for  Processor  IV  for  M  =  25; 

Known  Signal  Frequency 


PF 

X(PF,  M  =  25) 

10" 1 

1. 28597 

10-2 

2.  31209 

!0"3 

3.  03814 

io"4 

3. 61396 

Approximate  Detection  Probability 

The  RV  J  in  (C-19)  becomes  approximately  Gaussian  for  large  M  ,  and  we 
need  evaluate  only  its  mean  and  variance  to  approximate  its  PDF.  From  (C-24) 

and  (C-25) 

F'fo)  =  ,  E(»)  -  JV^(-^)t,M  £  .  <M3> 
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Therefore  (reference  6,  6.618  4) 


Eh)  '(j),>M  exf 


Also, 


(-my+)[^(^A)+x,(«y^  *  ">*  • 

•f"(o)  =  -t(5x)  =  j  Vir  exp(-  r")x,  (>»r)  (-i), 


yielding  (reference  9,  9.  6.  26  and  10.  2.  13  and  reference  6,  6.  618  4) 

E(A)  =  f  Jr  r  «p(-!di2l)[l„M-  ^rl_ 


--  \-iir  <*pKA) T*(MV+) 

*  |--^[''e*pW2^i 


The  variance  of  RV  5  is 


1  sc~ 1 

m5  s  <rs 


Vor(s)  =  I  -  t!?  ['~  e^ph^/h]' 

using  (C-36)  and  (C-34). 

Then  the  PDF  of  RV  X  in  (C-19)  is 

p(i)  -  Ae*r[-7PR; 


and 


?rot>(i>T)  a  lU*  f{  \W« 


(C-34) 


(C-35) 


(C-36) 


(C-37) 


(C-38) 


(C-39) 


Denoting  the  approximate  detection  probability  by  %  and  using  (C-30),  we 
have  v  \ 

where  Wls  and  q  are  given  by  (C-34)  and  (C-37).  respectively.  .Vs  the  sig¬ 
nal  strength  approaches  zero,  Pi^O  ,  and  (C-40)  approaches  the  approx, matt. 

false  alarm  probability,  ?F  :  f  \  (C-41) 

x  '  =  5(- v 


%  =  Uw=o>  =  I  - 
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using  (C-34)  and  (C-37).  This  result  is  identical  to  (032),  which  is  consistent 
with  the  observation  that  both  are  valid  for  large  M  .  Thus,  (040)  and  (041) 
constitute  the  large- M  approximation  for  detection  and  false  alarm  probabilities. 
For  small  M  ,  the  factor  WsAs  in  (040)  is  approximately  ;  comparison 

with  (05)  indicates  that  Processor  IV  is  1.  05  dB  poorer  than  Processor  I  for 
large  M  and  small  SNR,  w\. 


To  obtain  exact  receiver  operating  characteristics  for  this  processor,  (029) 
should  first  be  solved  for  the  required  X  =  X(PF)m).  Then,  that  value  of  \  , 
along  with  (030),  should  be  used  in  the  exact  relation  (026)  to  obtain  .  If, 
instead,  approximate  values  of  are  acceptable,  (040)  should  be  used,  where 
the  values  given  above  are  used  for  X(Pf>M).  That  is,  the  values  of  X  used  in 
(040)  should  not  be  obtained  from  the  approximation  (041),  but  from  the  exact 
relation  (029). 


C.  1.  5  PROCESSOR  V:  UNKNOWN  SIGNAL  PHASE; 

PHASE  SAMPLES 

From  table  1  and  equations  (5),  (015),  and  (020),  the  decision  variable  is 


JL  - 


1 

Pi  +  Qk  +  i  kK 

[(wi+-a^  +  \>l] 

E  U+IV 


(042) 


A  factor  of  l/M  ':as  been  supplied  for  convenience;  it  can  be  absorbed  by 
threshold  adjustment.  We  have  not  been  able  to  derive  a  tractable  method  for 
evaluating  the  detection  probability  exactly,  but  an  exact  false  alarm  probability 
calculation  is  possible  and  a  good  approximation  for  the  detection  probability 
has  been  attained. 

Exact  False  Alarm  Probability 

For  M  =  0,  we  have 

M 


i  =  iv| ,  >  v“  M-  ^  (a; Hi*)*  '  <c-43> 


The  second-order  CF  of  RVs  U.  and  v  is 


-{ff  dk  e*P[i 


2tt 


(C-44) 
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using  (C-18)  iund  changing  to  polar  coordinates.  Then,  the  second-order  PDF  of 
u  and  V  is 


Therefore 


pkv)  =5iy]Pf<)7  (C_45) 

=  ^  fj  &  r  T„(r  /uFT")  T*  (r/M  ) 

|-  PF=  Prot  [■$  *  t)  ~  if  J«dv  Jr-r 

0  (C-4 

-  i>  t,  w  j;  (*v) , 


{C— 46) 


upon  interchanging  integrals,  changing  to  polar  coordinates,  and  employing 
reference  9,  9.  1.  30,  Thus 


pF  =  i-  r* 


( C— 47 ) 


where  normalized  threshold 

X  3  V^T  (C-48) 

Since  Pp  depends  on  only  M  and  X  ,  Processor  V  is  a  CFAR  receiver.  The 
reason  for  the  threshold  scaling  in  (C-47)  and  (C-48)  is  that 


?F  -  I"  I0  dx  T,  (x)  exp(-  2^)  85  ^ **'  (c-49) 

or 

pF  -  enpl-xVa)  as  M  °4,  (C— so ) 

where  we  have  employed  (C-31)  and  reference  6,  0.618  1.  Therefore  the  choice 
of  X(pF;M)  in  the  exact  relation  (C-47)  becomes  relatively  independent  of  M  for 
large  M  .  Table  C-2 gives  values  of  X  ,  as  determined  from  (C-47),  for  M-2r. 

Table  C-2.  Thresholds  Required 
for  Processor  V  for  M  =  25; 

Known  Signal  Frequency 


PF 

X(PF,  M  25) 

lO"1 

2.  11272 

in-2 

2.99424 

10“ ;i 

3.  62192 

io-4 

4.  12829 

69 


TR  4529 


Approximate  Detection  Probability 

rv,P  RVs  u  and  V  in  (C-42)  tend  toward  Gaussian  RVs  for  largeM  .  We 
Fae  RVs  U  \  nf  u.  and  V  can  be  approximated  by  ajomt 

assume  that  the  second  or  ,,  «  t .  and  second-order  moments 

Gaussian  PDF.  Therefore  we  must  evaluate  all  first-  and  second 

of  u  and  V  .  From  (C-42)  we  have 

,  ,  ,  v=  r  I?; - ’  <c-51> 

u  =-  f/  a1  v  M  *=‘  J 

n  +  ()*)  +  Ph  J  L 

A  comparison  of  (C-51)  with  (C-21)  and  the  use  of  (C-34)  shows  that 

E(u)  ={ft  Me>op(-'H/,t)EPK/+)+T(B,,A)J  (C- 52) 

Also,  using  (C-18),  it  follows  that 


Also,  using  - -  . 

,  fr  i  ,1  b  I  ^  (C— 53) 

E(v)  =JJbO  cib-~ -  TfR  2m  *  2.  ) 

[_(W  F  «)  F  P  J 

since  the  integrand  is  odd  in  k>  . 

The  correlation  between  U  mid  Vis 

i  (yvi  V  Q«)  _ _  -  0  54  ^ 

eH=  ^$7,  Ell(mt  ,„>!.:?[(«+ «,)’+kT)  ‘ 

because  odd  integrands  cause  all  averages  on  RVs  ^  to  be  zero. 

The  variance  of  RV  v  is^  ^  ^  ^ 

=  Eiv‘)  ’  l  [(.«+ %)l+  IvT'[(w+%)' + k’TA  1  <t>55) 

-  m  Eli^Wl +  o-  . 

14.1,  w-t  t- farms  from  the  k*J  terms.  The  second 
tZZcZ  is.^iiiordiDg  to"(C-63).  -o;  and  the  first  term,  upon  changing 
to  polar  coordinates,  is  ~i 

AA  S'  1  A  ...  *  1  ! 


(C-55) 


coordinates,  is  — i 

r>  =  J-Jjrr(jei  * e  e*r[-i(',->'r«»+*Vj 

?•  .  ,  /  r»+..n  t  i- <txp(-"y4) 

"  y\"  J  ^  vvT  W*  m  V  2  J  H  Wj1 


(C-56) 
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where  we  have  used 


^  [do  s\n  0  exp(*  Ct>5&)  -  co3^) tXf(*  cc5B)  -  I ^(a)-J^(ci)~  (c_57) 

2ir  >r 

and  reference  9,  9.  6.  27,  9.  G.  2G,  and  10.  2.  13,  and  reference  G,  6.  G18  4. 
Finally,  the  mean-square  value  of  RV  u  is  given  by 

n  i  '5'rJ _ Qm^kY*vi+^)  1  ,  J  .jB±fl)a  ■  /._ 


c/,'  JL^E  _ _ L_LE  — U/|-i-Y'  -  -~u  (C-5 


The  first  average  in  (C-58)  may  be  expressed  as 


a-Ti-d  ^  T1  r  i  fi 
M  1  [  +  b1  3J  M  In'  J, 


(C-59) 


upon  using  (C-55)  and  (C-5G).  The  second  average,  by  inspection  of  ( C- 5 1 )  and 
(C-32),  is  obviously  W\*  .  Therefore  the  variance  of  RV  U  is 


(T1  =  JL_f.  _  ±exg(-y/i\  _  vyly  I 

°u  M  L1  v**  UJ  • 


(C-GO) 


It  is  immediately  obvious  from  (C-5G)  and  (C-GO)  that  the  variances  of  R 
and  V  are,  in  general,  unequal.  Thus  for  the  second-order  PDF,  we  have 


Then 


pM  S  (*’W '  . 

VB  =  ?™k  (i>r)  -  pa(AFc’>r)rj[  JuJvpP,') 


(C-Gl) 


U4v>T 


(C-62) 


U>/  >T 

Equation  (C-G2)  gives  the  approximate  detection  probability  ;  its  evalua¬ 
tion  is  in  terms  of  integrals  of  elliptieally  bivariate  Gaussian  functions  over  off¬ 
set  circles,  for  which  there  are  few  tables  available  (e.g. ,  see  reference  10). 
Accordingly,  a  different  approach  is  used  here;  i  e. ,  the  integral  on  V  in(C-G2) 
is  first  evaluated  using  (C-4).  Then  there  follows 

st  c’A  r  .71  -r  /  /  m, 


g  =  dreip[_+(X-t)r|[2|(i/rrA‘')-l]i 


(C-G3) 
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where 


,  _  ty,  n  _  _xl  _  d  =  — 3 

®  C  ;  ^./a  ~  2Mflvl  '  °V* 


(C-64) 


Equation  (C-48)  has  been  used  in  (C-64);  the  quantities  required  in  (C-64)  are 
given  in  (C-52),  (C— 56),  and  (C-00).  Numerical  evaluation  is  easily  effected  by 
using  approximations  to  ^  (e.  g. ,  see  reference  9,  26.  2. 17). 


As  the  signal  strength  tends  toward  zero,  there  follows  from  (C-20), 
(C-52),  (C-56),  and  (C-60), 

<Tu  —  C2M)'- 

Therefore  an  approximation  to  PF  is  afforded  by  (C-63): 

PF  =  PD(^-0)  =  <E*p{-  \7^), 

which  is  consistent  with  (C-50)  for  large  M  . 


(C-65) 


(C-66) 


Equations  (C-47)  and  (C-63)  constitute  the  analytical  results  for  Processor 
V.  Values  of  \  are  determined  from  (C-47),  and  then  substituted  in  (C-63)  and 
(C-64). 

A/ 

It  is  of  interest  to  evaluate  for  small  SNR,  and  to  compare  it  with 
Processor  II,  which  also  uses  amplitude  samples.  From  (C-52),  (C-56),  and 
(C-60) 

wv'fc'fv  U  *-*•>.  (C_67) 


Therefore  'C-62)  yields 


yieiufa  _ 

%  ~  ij  &  dv  W  2M  cxpf  m(u-Vt  ^  -Mvj 

~  Q  m,  /IFt)  =  C^Vm1  ",  x)  Os  w  • -*  ‘ 


(C-68) 


where  we  have  changed  to  polar  coordinates  and  ased  (C-8)  and  (C-48).  A  com¬ 
parison  of  (C-68)  with  (C-9)  reveals  that  Processor  II  functions  equallv  as  well 
as  Processor  V,  while  using  a  value  for  )vt  that  is  \/?/2  less.  Thus,  Processor 
V  is 

20  log  (^")  ~  l-OS"  dB  (C-69) 

poorer  than  Processor  II  at  low  SNR  and  large  M  ,  where  the  central  limit 
theorem  is  valid  for  the  analysis  of  Processor  V.  Thus,  loss  of  amplitude  in¬ 
formation  does  not  cause  much  degradation  in  detectability;  phase  information 
is  the  major  factor  in  determining  detectability.  It  will  be  noticed  that  the  amount 


- -i-w-'n-mniy jrtiii ii ■ ,  -.ijr;,; _ 
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tiie  performance  of  Processor  V  is  below  that  of  II  is  exactly  the  same  as  that 
by  which  Processor  IV  is  below  I.  These  are  analogous  processors,  in  that  each 
pairing  differs  only  in  its  use  of  amplitude  information;  they  make  the  same  use 
of  phase  information. 


C.  1.  6  PROCESSOR  VI;  UNKNOWN  SIGNAL  PHAS1  ; 

FITTED  PHASE  SAMPLES 

The  decision  variable  for  this  processor  (table  1)  was  given  as 


M  K, ,  K  V  M  K-i  ' 


(C— 70) 


We  are  unable  to  exactly  evaluate  or  P  in  ;uiy  simple  form.  Rather,  we 
outline  a  numerical  approach  that  ould  be  used, and  then  derive  approximations 
..0  the  detection  and  false  alarm  probabilities. 

Exact  Detection  Probability 


The  CF  of  R V  X  in  (C-VU)  u 


p^)exp  ^ a*]  exP  [~l  In9*)  _  . 


(C-71) 


Now  the  second  exponential  in  (C-71)  can.  be  simplified  by  means  of  the  following 
artifice  (reference  11,  (22)): 

exp(-<ryy2)  =  f dv  (yn-T-1)  exp(-jp-  +  '/*x,i.  (C-72) 

The  quantity  p’  on  the  left  of  (C-72)  is  replaced  by  ix  on  the  right.  Identifying 
'  .  M  r — \ 

y  =  ^  M  =  Vl'2  -  l+l,  (C-72) 

the  multiple  integral  for  the  CF  in  (C-71)  c;m  then  be  expressed  as 

e(„\  _  I  A. _ WaNSV  L  ^  —  T -1-2-^  I  Jj/Lvril^ 


(C-72) 


-f(f)  -  dx  wpOxY^fda  pfe)e*pl/if^  +  11 


I 7 

iMLL  expfe’  JjiJe  ffs)  e]( 

l+l  J  V  c  /  (J  L  J3 


(C-74) 


TWW  ^^rw5»',iujp)i,!^i  w  >|  lOfnacMdUP’ 
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For  a  particular  v;ilue  of  %  ,  either  form  of  the  inner  integrals  on  £  in  (C-74) 
can  be  quickly  evaluated  for  many  values  of  x  by  means  of  an  FFT,  permitting 
immediate  evaluation  of  the  outer  integral  on  X  .  The  PDF  to  use  for  ylo)  is 
given  in  (A- 19)  and  depends  on  ^  .  rFhe  final  evaluation  of  the  cumulative  dis¬ 
tribution  of  RV  X  i3  available  by  using  the  methods  in  references  7  and  8.  The 
exact  Yr  may  be  evaluated  in  a  slightly  simpler  fashion  since  pfG)=^Tf^  ^ 
in  this  case,  and  the  integral  on  Q  in  (C-74)  is  the  error  function  of  a  complex 
argument  (reference  9,  chapter  7). 

Approximate  Detection  Probability 

The  RV  X  in  (C-70)  i  proximately  Gaussian  for  large  M  .  Therefore 
only  its  mean  and  vari;mc(  must  be  computed  to  obtain  approximations  to  7^ 
and  PF  .  Since  PD  depends  on  the  precise  value  of  1^0  ,  we  will  select  the  best 
value  to  maximize  For  the  test  in  (C-70)  where  |6k|  x  tt,  the  best  value  of 
(j^  is  zero.  This  point  is  discussed  further  in  the  main  text. 

in  order  to  evaluate  the  mean  and  variance  of  RV  X  in  (C-70),  we  will  need 
the  moments  of  RV  9.  For  small  SNR,  we  use  the  approximation  given  in  (A-20) 
to  order  T^/tr3-  in3  ;  we  obtair. 

^  =  E(en)  =  fefe  e*  p(&) 

1  o  /  /  o  n  r  \ 


(  de  6  ’  [| ^  Co  s  e  -V  cos(2e$|  t  io  Order  n\2 


(C— 75) 


It  follows  that  /A,  =  o  for  vi  odd,  which  is  the  result  of  our  choice  of  ^  =  0 . 
Using  integration  by  parts  for  A  even,  we  obtain 

,K>  = 1  > 

_  -3|c  ftp  4  \  h)3;  Xt  orefer  Y*  ,  (C-76) 

^  -  2X'frr'(vX~t)'rt  4-  ^  (2^-3)  ™7;  h>  order  m2- 

The  mean  of  RV  X  is  obtained  from  (C-70)  as  follows: 

Eli]  -  =  rX- m)  -  •  <c~77’ 


The  square  of  RV  X  can  be  expressed  as 


M  rj_ 

Jt^  -  rj  +  jf  2>_J^G»0r,©F 

*  M  ^Ti  K  <V  w-i  n  k#m<N/P-i 


(C-78) 


(C-79) 


The  means  of  each  of  the  three  summations  in  (C- 78)  are,  respectively 

N i/v+ 


The  mean-square  value  is  then  available  from  (078)  and  (079),  and  the  variance 
of  jP  is,  finally, 


Vnrfi}  =  e£T]  -  e"U]  =  2ri"  mJ [v  A  £t]  =  •  'C'80* 


The  ratio 


^‘-VS’-F 


5td.  Pev.  ^  e1} 


■fov-  lam?e  N. 


(081A) 


(081B) 


The  last  relation  would  be  obtained  from  (070)  by  simply  dropping  the  last  term. 
Thus  an  approximation  to  Processor  VI  v/ould  be  afforded  by  ~jvf  ^  6\<  ;  this 
possibility  has  not  been  pursued  further,  however. 

The  probability  of  RV  Jl  not  exceeding  a  threshold  T  is  approximated  ac¬ 
cording  to 

?rok(i  <  t)  s  («  (hrQ1)lt«r[-  =  f  pip)  -  ^ 

—  cG 

by  using  (04).  1’he  quantity  is  given  by  (C-81A),  (077),  and  (080)  in 

terms  of  p-x  and  j. ,  which,  in  turn,  are  given  by  (076),  or,  more  generally, 
by  (075).  For  large  M  (and  -6),  we  obtain  the  approximations 
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where  A  E  ~^r~  (j  -  :~*  1")  .  A  comparison  of  (C-83)  widi  (C-5)  shows  that,  when 

vll  =0  ,  Processor  VI  is 
To 

20  lo^  (l/ ?^i)  -  1-3*1  <lB  (C-84) 

poorer  than  the  optimum  phase-coherent  processor  (I)  for  large  M  and  small 
SNR,  whereas  a  comparison  of  (C-83)  with  (C-4b)  et  seq.  shows  that  Proc¬ 
essor  VI  is  0.  34  dB  (1.  39  -  1.  05)  poorer  than  Processor  IV  for  large  M  and 
small  SNR.  Thus,  for  known  signal  phase,  ,  the  nonoptimum  combination  of 
phase  samples,  (C-70),  is  only  0.  34  dB  poorer  than  the  optimum  combination 
(see  table  1)  for  large  M  and  small  SNR.  However,  these  comparisons  may  not 
be  relevant  for  the  range  of  useful  and  T£>  ;  this  relevance  can  be  determined 
only  from  the  detailed  receiver  operating  characteristics  presented  in  the  main 
text. 


Ek 


C.  1.  7  PROCESSOR  VII:  UNKNOWN  SIGNAL  PHASE; 

PHASE-DIFFERENCE  SAMPLES 

The  decision  variable  for  this  processor  (table  1)  was  given  as 

Re  £  evpO'O*- /#„_«)].  (C-85) 

We  actually  consider  a  more  general  processor  here,  namely,  one  that  uses 
weighted  phase  differences 

Jt  =  (C-86) 

because  better  performance  may  be  attainable  for  jxto  man  for  ^~o  .  In  trade, 
however,  the  CFAR  capability  is  lost  for  ptQ. 

Exact  Detection  Probability 

The  one  case  in  which  a  closed-form  expression  for  the  OF  of  RV,I  in 
(C-8i3)  is  attainable  is  where  yx-l  .  To  show  this,  we  use  (5)  and  (C-15)  to  ex¬ 


press 


^  K-ij  =  Jh  +  A  . 


(C-87) 


The  CF  of  Jl  is  given  by  the  product  of  the  CFs  for  IlVs  and  since  fa#} 
and  {\>^j  are  independent  (see  (C- 18)).  Also,  the  CF  for  jt^  is  a  special  case 
of  that  for  JA  ,  obtained  by  setting  vn  =  0  in  the  quantity  Ja  .  In  order  to  obtain 
the  CF  for  RV  Aa  ,  we  use  reference  .1,  appendix  C:  Matrix  P  has  elements 


£ 


= }  rt,  r,-j!  = 

0j  otherwise 


I  s  i,j  *  M, 


(C-88) 
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and  all  thr  elements  of  eolumn  matrix  M  are  equal  to  M  .  The  covariance  ol' 
RVs  is  equal  to  the  identity  matrix;  therefore  matrix  ,A  is  equal  to_P 

(reference  1,  (C-H)).  Let  the  normalized  modal  matrix  and  characteristic  value 
matrix  of  B  be  denoted  by  Q,  mid  X,  respectively,  Then  the  CF  of  Ja  is 


(C-«9) 


where  matrix 


^  ;  (S  w  =  ynfi  1  , 

It  follows  that  the  CF  of  KV  J \,  is  given  by 

ft1- 

mid,  therefore,  the  CF  of  RV  JI  in  (C-S7)  is 


Evaluation  of  the  characteristic  values  [ \k"j  of  Jg.  is  possible  from  ref¬ 
erence  12,  ((57)  and  ((is),  and  are  given  by 

\  =  1  £ks  M-  tc‘i,:)) 
KvaL  ation  of  the  normalized  modal  matrix  0  of  matrix  .g  in  ( C—  88)  requires 
cor  i|  ater  computation,  and  (C-92)  c;m  be  evaluated  only  with  such  aid.  The 
cumi  lative  probability  distribution  and,  hence,  the  operating  characteristics 
of  Processor  VII  for  ^.  =  1  in  (C-8(i)  are  then  available  by  means  of  the  techni¬ 
que1:  described  in  references  7  and  H.  This  exact  approach  to  the  detection  prob¬ 
ability  for  I  has  not  been  pursued  further  in  this  report. 


(C-90) 


(C-i)l) 


(C—  DC) 


Exact  False  Alarm  Probabilities 

fo r  -  I  .  The  false  alarm  probability  for 

able  from  the  above  results.  For  )vi=0,  (C-90)  yiel 

M 


The  false  alarm  probability  for  u  =  |  is  immediately  avail- 
/e  results.  For  )>i=0,  (C-90)  yiekls  f*K-C,  :md  (C-92)  becomes 

M  vt 

fir)-  7T(l-i2V?)  .  (C-94) 


I'he  PDF  of  RV  Ji  in  (C-S7)  is  th(-n  given  by 

frVj]  q  /  j 

pW '  i  X  "PC  fS 


(C-94) 


(C— 913) 
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where 


M 

c„  =  TT 

K  -  I 


X 


m. 


<  m  <  [m/X]  ( 


(C-96) 


it=,  aw"^k 

;Uld  {yj j  are  gWen  in  (093).  The  probability  of  false  alarm  is 

PF  (pi)  =  TWlC^r)  =  jT<W  ty) 

fc,(pa)  WT,a 

V*  .- 1 


(C-97) 


it  is  convenient  to  define  a  normalized  threshold  as 


in  which  case  (C-97)  becomes 


er  » 


(098) 


iv\  -  i 


'_  4CW  x 

co5(mw) 


( C— 99) 


The  reason  lor  the  scaling  in  ,o9»,  is  that  Pf  in  JC-W  “* 

as  M-.  TO  son  this,  .0  can  either  use  (C-»7,  and  note  that  a  appioaches^a 

Gaussian  RV  with  a  zero  mean  and  a  variance  ol  <T  4^')  or  B 
can  investigate  (C-94)  for  large  M:  We  have 


f(?)  ' 


(C-100) 


_ ___ _ I _  -for  5h»va ))  !x, 

,7)775  A„  +  (■  lit)’  Y¥/>  y  ^ 

But  the  sum  of  chararteristle  values  us  aero  because  the  diagonal  elements  ol 
V  in  (088)  are  all  zero.  Also  /  m  \i  O  , 

2.\i=  (-2\)  -  Jo 

^  ^  *=•  V*-*  /  *•'  (0101) 

K~  i 


using  (093).  Therefore 

f(T) 


r  4v  \0X1  r  a «  (Oio2) 

_ =  e*p  -<r\M-l)j  ]  Tjr  5«ta/)  T, 

iL/  v 


fdr  i  -wo-^-O t;‘  1  u  v  n  u  <• 

u  rs  •  n  rr  With  I  zero  mean  and  a  variance  of  2<rfyl-i).  Ihcrf>forc 
;f;  ii:  «.«.«  becomes  relatively  independent 

Ir  m!  tor  large  M  Table  C-:«  gives  values  ol  X  ,  as  determined  Iron,  (C-..J). 

or  M  if  and  ^  “ 1  • 


?S 


■ttVa .  .. 


rtf/i  -rilla'/ri1—  ■*■•£  iMli  tiiit  ,  v  c„,  ... 
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Table  C-3.  Thresholds  Required 
for  Processor  VII  for  M  25  and 
m-  1;  Known  Signal  Frequency 


Pp  (M=  1) 

A  (PF,  M  =  25) 

-1 

10 

1.  25712 

-2 

10 

2.  40560 

10“’* 

3.  35438 

-4 

10 

4. 21662 

Pp  for  p*=0  .  The  false  alarm  probability  for  j^=C  in  (C-86)  can  a 
derived  in  a  tractable  furm.  This  processor  is  the  phase-difference  pr< 


(C-H5),  that  we  express  as 

M  \ 

(^co5  Q,  -05  0K.(  +  Sno  Sin 

The  CF  of  Jl  is  then 

-fff  )-  £:jexp[if  ( 105^0)5  8^,  +  5|/)0*5M0K.,)]j  , 

Now  the  average  on  Q(  in  (C-104)  is  (see  (C-15'.  through  ( C  - 1 R ) ) 

— _  {  <lG,  expjjf>(cc6  0ico5  01+-5in0J.Jtio  0,)  =  T0(j)  , 


also  be 
rocessor, 

(C-103) 

(C-104) 

(C-105) 


A  similar  procedure  applied,  in  turn,  to  0j  0  0  yields  the  closed-form 

CFfori:  ’ 

Iff)'-  (010(5) 

Then,  according  to  reference  7,  (7),  the  exact  false  alarm  probability  is 
f>b  =  o)  -  ?rol,^>T)=  i--f  iiffT) 


r(^ 


~S— »  JL 

K-l  \ 


(C-1071 


wh'»re  no’mali/.ed  threshold 

X  = 


(C-108) 
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By  an  argument  similar  to  that  in  (C-28)  through  (C-32),  we  find 

PFl;*  =  o)~  cF(-\)  **  (C-109) 

Table  C-4  gives  values  of  X  ,  as  determined  from  (C-107),  for  M=25~  and 


Table  C-4.  Thresholds  Required 
for  Processor  VII  for  M  25  and 
M  0;  Kno'.vn  Signal  Frequency 


n 

I5F  (h  -  0) 

X(PF,  M  =  25) 

-1 

10 

1.  28615 

10"2 

2.  31149 

10~'5 

3. 03591 

10-1 

3. 60941 

Approximate  Detection  Probe  hility 


The  RV X  in  (C-80)  for  large  M  is  approximately  Gaussian.  Accordingly, 
we  need  evaluate  only  its  moan  and  variance.  By  using  (C-15)  and  (C-20),  we 
first  express  (C-86)  as 


11 

(vv\  -p  M  \\  ^X-l 

x-  [(»+««)’+  mTLmj**  aj,T 

(C-110) 

where 

v  - 

(C-lll) 

Then  we  define 

c,  * 

M+  Qk  S  - 

(C-112) 

[(ui+QK)2  +  KY  '  K  [CwtFflJ  T  ^kT 

and  express 

Jl  =  C 

(C-113) 

where 

U  "  M-l  ^  (CkQ-i  + 

(C  114) 

HO 


-  -  Sisiu 


\ .. 


Til  -152!) 


At  this  point,  we  derive  several  generic  results  for  the  ItVs  ^Q,j  and  [S*} 
that  we  will  need  both  now  and  later.  (The  derivation  of  the  mean  and  variance  of 
U  resumes  at  equation  (C-121).  )  We  shall  use  (C-1H)  et  seq.  freely.  We  have 
the  mean 


c* 1  4  (iJXW  £  e*r  (' 


:  r0^kT  ^  n  2  > 

.  fir  V-  stride  r’^cose  «J-  xlilapW 

0  2ir  L 

=  f° dr  v1  2%xj>(-  J— • X,  (wr) 

=  2  r(2-v,)‘v'  exp(-wiV2)  (Ff  (2-V;  2 ;  ^  I2-) 


(C-115) 


pzl 

**  r, 


(3_^_)  exp^V^.f;  (^;  2;  £)  *  C , 


where  we  have  transformed  to  polar  coordinates  according  to  «=  r«>s0- ta,  t = v~si»^f 
and  used  reference  0,  8.431  5,  6.  031  1,  and  (C-lll).  Also, 


=  j  cjr  r  de  r  exp[-  "X? -  =  0,  (C-iiG) 

rr\  ~>-rr~  ~ 


using  the  oddness  of  the  0  -integrand. 


The  cross-moment 


%50  Sin8ex£  J±-22n^056  (C-117) 


using  the  oddness  of  the  0  -integrand.  Also, 


Sk  =  i0  dr  r  {c!o  r'^snd  exp[- 


-2^yircosO-fW 


f°,  3-“V  .  /  r1+hi1\  ^(wr) 

-J0^r 

x  P  (2-2»;  e*P  (-W.V2)  ,ff  (2  -  2P ;  2 ;  vX/z) 

-  2^  r  (i+p)  ^(-wV2)  0  V ; 2;  E  ; 


(C-118) 


TrffireegWBPH 


P  '• 


4 


'4,  *■ 
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using  reference  6,  8.431  3  and  6.  631  1,  and  (C-lll).  Finally, 

_  $  f  ..  2-4y  J,  r  \r’x--2.n\rc(&6+vn  I  fC-1191 

C‘  =  j<!rr  r  C05  6  «*[_-  ^  J.  ' 

7-  |  _■  u  ;n  ,0119)  and  use  reference  6,  6.  631  1,  and  (0118) 

We  express  Cos  6  =  |-si*  0  m  (L  UJ)  a,lu  uo^ 

and  (C-lll)  to  obtain 


r;-5>  -^«p(-^)t0M-^ 

...  r  T-V  XT.,  ;  /  r1 


have 


(C-120) 

Now,  we  resume  the  calculation  of  the  statistics  of  RV  u  in  (C-114).  We 

(C-121) 


_ i  _ * 

lx  -  C  +  5 


C  =  ™u, 

using  (C-115);  the  second  moment 


-110)5  Ulw  bucuuu 

—a  -  -L_  ^eCckCk.,  +  5,5k.,)(9  <*.,+  ^  5--')] 

W  "  (M-0 


(fi-ij 


4~ 


4- 


±e(CC,.,  +  5,5b.,T 
Z  e((c„Cb.,  +  ^  5,-')(c^  c" + 

<-'i 

j  (cK  Q., + 3,  S*.,)  ( C,  Q->  +  3,,  WOj 
Ej(Ci<CK-,  +  5sSK-Xc<.Ci+  ^ 


M 


(C- 122) 


K=3 

H 


4  „ 

K,JP-2 

The  four  averages  in  (0122)  are  given,  ivspectivelyyliv 

cj  Qt  +  2  S„  C„  SK„  C,-,  +  c  +  5  .  , 

c>-  Z,  C=  +  ^  ^  ^ 5-  ^ +  s'2"_ _4  .  _.,c  123, 

-.-rr^fT  5Z+5,sC^,C1^5:,5<-,-c  c, 

CK  c-K-(  CMb44-)-V,^ 

_A- 

c, 


82 


upon  using  (C-115)  through  (C-120).  Combining  (C-121)  through  (C-123),  it  fol¬ 
lows  for  the  variance  of  U  that 


O'1  -  — 
lA  M-l 


"(?  +  c/  +  sMc-  i£rC  (c1-  0], 

The  PDF  of  U.  is  then  approximately 

1  (u-y%)*~] 


2<rJ 


and 


P(u)  ? 

PD  =  P«S>P>t)  =  R-»l(u> 
~  ([) 


Mu  T 


(C-121) 


(C-125) 


((’-12(5) 


The  quantities  W)M  and  are  available  from  (C-121)  and  (C-124),  upon  em¬ 
ployment  of  (C-115),  (C-11H),  and  (C-120).  The  approximate  false  alarm  prob¬ 
ability  is  obtained  by  setting  in  the  above  results  and  is 

%  =  £(.->>),  <c-1271 

where  normalized  threshold 

^  2.r'  V[\+y)<rV  '  (C-12,s) 

(This  agrees  with  (C-9H)  and  (C-108)  for  -  !  and  J.<  - 0 ,  respectively.  )  'Flier. 
(C-120)  becomes 


*  t-  'r(M  \\ 

-x;. 


(C-129) 


Equations  (C-127)  and  (C-120)  constitute  the  approximations  to  ft  and  FJj  for 
arbitrary  u  ■  Threshold  X  is  determined  from  (C-127)  for  a  specified  fp  , 
and  then  employed  in  (C-120)  to  evaluate  P0  .  However,  for  jA  =0  ov  J4- 1  ,\ 

should  be  determined  from  (C-107)  or  (C-00),  respec  ,ivel\ ,  and  substituted  in 
(C-129)  for  evaluation  of  ?D  .  (Of  course,  the  best  approach  for  -  I  is  to 
evaluate  the  exact  Fp  via  (C-92).  )  The  quantities  Mu  and  cr^  are  available  in 
(C-121)  and  (C-121). 

Approximations  for  Small  w 

For  small  SNRs  per  sample,  Vt\«  I  .  In  this  case,  the  following  approxima¬ 
tions  result  to  order  w  : 
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C  ~  2^P(^)«, 
s3-  2r'r(i+/H), 
C~2r'  p(i+r), 

m^V^FV, 


(0130) 


Then  (C-129)  becomes 


pM3^) 

rO+y) 


yy\~  -  >s  w- 


(C-131) 


Two  important  features  should  be  noted  from  (C-131).  First,  the  dependence 
on  ya  is  according  to  wi1  ,  not  according  to  ft;  as  for  Processor  I,  (C— 5);  Proc¬ 
essor  IV,  (040);  and  Processor  VI,  (083).  This  result  is  also  discernible  in 
the  exact  CF  calculation  1092)  for  ^=1  ,  where  the  dependence  of  ^  is  accord¬ 
ing  to  vy\*  (see  (090)).  Thus  the  phase-difference  processors  have  a  small-sig- 
nal  suppression  effect  for  all  values  of  weighting  constant  /a  in  (086).  Second, 
tuere  is  an  optimum  value  of  yu  to  maximize  the  coefficient  of  hi1  in  (C-131). 
The  calculations  shown  in  table  05  indicate  that  the  best  value  of  weighting /t 
to  use  in  (080)  is  unity;  however,  the  factor  has  a  broad  maximum  about  the 
pointy  =  !  .  (Also,  the  best  value  of  p  for  larger  SNRs  may  not  be  1.  )  The  loss 
i.i  detectability  at  JJ-- O  versus  p-\  is  available  from  (C-131)  and  table  C-5  and 
is  10  log  (4/rr)  =  1.  05  dI3  for  small  SNR;  once  again,  the  loss  cf  amplitude  infor¬ 
mation  incurs  a  1.  05-dB  degradation  in  performance  for  small  SNR  and  large  M  . 


Table  C-5.  Multiplying  Factor  Used  in  (C-131) 

—  n  :  r“m . r 


1/21  1-1 


/r(f  '  m)|  *74  .  785  I  .  953  ,  1 


566  .884 


iViiSiv*  j.Vi.fci- 


-  i>. .  .i  tw-vii li  W 


C.  2  UNKNOWN  SIGNAL  FREQUENCY 


In  accordance  with  the  discussion  at  the  end  of  section  2,  the  only  processor 
we  will  consider  here  is: 


C.  2.  7  PROCESSOR  VII;  UNKNOWN  SIGNAL  PHASE; 

P  H A  S  E -  DI F  F  E  R  E  N  C  E  S A  M  P  LE  S 

The  decision  variable  for  this  processor  (table  2)  was  given  by 

i  M  ,i 


i-  ^  e*p(ie-,<-  'O  . 


1C-132) 


We  actually  consider  a  more  general  processor  here,  namely,  weighted  phase- 
diffe  rences 


J  -  ^  e*p  (<  '  ^-i) 

K-2  ) 


VC-133) 


because  better  performance  may  be  attainable  for  than  for  /<.  =  0  .  An  exact 

relation  for  the  detection  probability  is  not  available  for  any  value  of  u .  (For 
\  ,  the  fundamental  problem  is  evaluation  of  the  PDF  of  »  where 

are  complex  independent  Gaussian  RVs  with  nonzero  means.  )  We  will, 
however,  obtain  approximations  to  the  detection  and  false  alarm  probabilities 
for  general  values  of  pt  . 

Exact  Fals..  Alarm  Probabilities 

In  a  manner  similar  to  that  given  in  (C-110)  through  (C-114),  we  expi’ess 
(C-133)  as 

X  -  <r  }  (M-i)I  UAi  v|  =  o  J  M  4- 1  v|,  (C-134) 


where 


,M-i  k 


g-(C,  C*,  +  Sk.,), 

k  (s*c„  -  c« 


(C- 135) 


V  -  Ml  X  X-K-l 
The  second-order  CF  of  RVs  u  and  V  is 


=  Ej 


ttf,'})  =  E[e^p((?H+iyY)] 


(C-13G) 


TR  4529 


Now  if  we  define 


c'  •-  (?C,  +  7  5.)/(«-') ,  P  ={?S1 - 7  G)/(m-i), 


the  average  on  C,  and  5,  in  (C-lIKi)  becomes 


E  £ Bvp  (  I  ex'  C  i  +  5,  )j 

=  Jfda  A  e*P[i  e*f  (-  ^T^-) 

-  |^dr  r  jcle  evp[i  (°f  to50  4-^3 >in  {>)]  ^P^-  r/2) 


Sr 


-  Or  r 

Jo  ' 

--  fdr  r  exp(-  r2A)  J0  (}^r 


(C-137) 


(C-138) 


where  we  have  employed  (C-tll)  and  (C-112),  set  m=0  ,  and  have  changed  to 
polar  coordinates.  At  this  point,  we  are  stymied  in  attempting  to  evaluate  the 
average  on  Cz  and  Sx  in  (C-l.'W),  except  in  the  two  special  cases  of  w-o  and 

rl- 

Pf  for  -0.  For  pt=o,  we  observe  from  (C-lll)  mid  (C-112)  that  QV!\=| 
in  which  case  (C-138)  becomes 


70{-^ir^, 


(C-139) 


A  similar  step-by-step  procedure  for  the  remaining  averages  in  (C-13G)  yields 


m-I  / 


fM=T 

Then  by  the  procedure  given  in  (C-  1  1)  through  (C- !(>),  we  obtain 

Pn,\>(i>T)  -PF(f  =  °)  =  I-  Edx3;(y)T0M^(x/T) 

=  I-  fdv  T,(x)  oT'CVsE i), 


(C-140) 


(C-141) 
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where  normalized  threshold 


\«VffrT  (C_1'I2) 

The  scaling  in  (C-142)  is  chosen  such  that 

?F  [jX  ~  o)  ~  e  Xp(-  X1/^)  Q5  M  -*  ■*>,  (C-14 

wliich  follows  when  (C-47)  through  (C-50)  are  used.  Table  C-G  gives  values  of 
X  ,  as  determined  from  (C-Tll),  for  M-25"  and  p  =0  . 

Table  C-6.  Thresholds  Required 
for  Processor  VTI  for  M  =  25  and 
M  0;  Unknown  Signal  Frequency 


(C-143) 


.  (M  0) 

X  (Pp  M  -  25) 

io"1 

2.  14258 

1 — 1 

O 

to 

2.  99250 

10-3 

3. G1779 

1— * 
O 

1 

4. 12110 

?F  for  n  =  I  .  For  I  ,  we  observe  from  (C-lll)  and  (C-112)  that  CK  -  4* 
and  5k  =  i)H  .  Then  (C-13  1)  and  (C-135)  become,  respectively, 

u  =  «2-(QkQk-i  + 


v  =  S  (W  ah  bK.,). 

K 

The  second-order  CF  of  U  and  V  is 

E  ^  e*f>  [i  f  X (Qk  4K-i  +  k  bK-i)  +  "/  kTI (bn  «k- i  “  bn-.)]] . 


(0144) 


(C-14  5) 


The  average  on  flj  in  (C-14 5)  is 


Jda,  e*p£(  *yf  bz)^(2v)  &.ip{-"/2)  *  b») ])  (C  14 

vhile  the  average  on  b,  yields 
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(C-147) 


(C-148) 


exp[-i(^-7aJ]. 

The  product  of  these  averages  is 

exp[-i(f2+v!)(^+^3 

(This  result  is  also  obtainable  from  (C-138)  for  =  I  ,  using  reference  6,  G.  G31  4, 
when  we  observe  the  scale  factor  f-H  used  in  (0134).  )  Now  the  averages  on  4, 
ar.d  bj  in  (C-145)  are  very  similar  to  those  in  (C-14G)  and  (0147),  except  for 
the  added  exponential  (C-148).  Generally,  the  average  on  3*  takes  the  form 
,  r  r  .»r  x  X,JfW)TI  ,  ,  >\1"A 

^V*P 


)oneniiai  uoiciaiij,  ^  ~  ~  &  t 

:  ===r<». p[-l(t<V.,+7U  xK(?V)]  ’  1  “  M  '' 


(C-149) 


where  TK{^)  is  a  polynomial  in  j  that  satisfies  the  recurrence  relation 

Xb)-  '■  .e-iso) 

In  fact, 

t  (^ .  ^  (V)  y !  =  T  ^  - M '  -  •  <c-151> 


K- 0 

The  first  few  polynomials  are* 

T0  k)  -  * 
t  (a)  =  i 
TG)= '  +a 
T(y)=i  +  2j 

ry(^= 


(0152) 


The  product  of  the  averages  on  QK  and  ^  in  (c_145)  1S  then  glVC”  by 


♦The  factorization  of  TM(y)  is  presented  in  (C-1G2)  and  (C-1G1). 
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_ I _ „T  jO*  1  +b  a)(fW)  1  ,  | s  K  s  M 

W+f)  ^  -r.if+tf  J  >' 

Finally,  the  averages  on  and  bM  in  (C-145)  are  given  by 

jL  jjX  dlM  evpf-K^U’Xif^i^]] 


feVm  T,_,(?vy> 


(C-153) 


(C-15‘1) 


_TM  (t+f) 

Equation  (C-154)  is  the  second-order  CF  of  the  RVs  u  and  V  in  (C-144).  Then 
by  the  procedure  given  in  (C- 44)  through  (C-4G),  we  obtain 

ProK^T)  =  PF  (/<  *  0  *  Pi-oj,(|ii-t-ip|>T/<e) 

,  l-_f  d*  a;(x)/TM(^xi) 


■  i  r°dif ,  (x)/ "Pd  (  2(m-i>xi  )  ■ 


where  normalized  threshold 


^  _  rzr. — o  * 


^2{ r 

The  reason  for  the  scaling  is  that 

PF  (^  = »)  ~  as  M 

To  prove  this,  wc  use  (C-151)  to  show  that 

T  /  X  '  .  ^  jMJfH  _  a-/-*!''* 


1  M 


2  ( M  -  i)  X1 


K-o 


-  ^  (M-2*)!(M-lT  *! 


(C-155) 


(C-15G) 

(C-157) 


(C-158) 


Therefore  (C-155)  yields 


PF(pl)~  \-f^x  ^Wexp(-f^)  =  exp(-xy2)  (C-i59) 


upon  using  reference  G,  G.  G31  5  and  8.  352  1. 
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Now  we  will  evaluate  (0155)  in  closed  form  by  ilrst  determining  die  roots 
of  the  poI\  nomial  TM(y)  .  We  notice  that  RV  W  in  (0144)  is  equal  to  RV  d/tP 
in  (C— 87)  for  W--0  ;  but,  the  CF  of  the  RV  jQr  is  given  by  (094)  and  (093). 
Therefore  the  CF  of  RV  u  in  (0144)  is  equa'  to 


M 


X-\ 


where 


W  v_, 

-  IT  (|++AV)  , 

K-i  ' 

(C-1G0) 

isK'M. 

_M4-I  /  9 

(0161) 

However,  from  (0154),  the  CF  of  RV  ft  is  given  by  Tpj  (tfi)  :  Set  y=0  in  (0145) 
and  (0154).  Therefore 

WtI 

T„(y)  =  TT  (i+ yi a).  <c-162> 

In  order  to  evaluate  (0155),  we  first  expand  Tj^'(0  in  partial  fractions: 

I  NO  « 

T“’(u)  =  -raji - - -  -  21  - TT-  ’  (C-163) 

a  if  (i+yla)  K"  i+Yr!* 

(Notice  that  -  1  since  Ifjo)  =  I  .  )  The  coefficients  in  (C-163)  arc 

given  by  K= 1 


[mM 


P 


I hi*) 


<^0S  ^ 

r  kT 

J 

(C-164) 


Kfd  ltd 

Substituting  (C-163)  in  (0155),  vc  obtain 

f' 


=  1 


ng  (C-163)  in  (O-LOU),  ve  ornain  ^ 

(/>=')  =  1  -  ^ [' + 
r  a! 


(C-1C5) 


w  here 


#-  l)/z 

A 


-X  = 


(C-166) 
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Hie  integral  in  (C-165)  was  evaluated  by  means  of  reference  13,  part  II,  521.  6. 
l’he  function  K,  in  (C-i05)  is  a  modified  Bessel  function  of  the  second  kind  of 
order  one.  liquation  (C  105)  is  a  closed-form  expression  of  the  false  alarm 
probability  for  ^u-l  ;  the  quantities  [”RK^  arc  given  by  (C-104),  and  the  quan¬ 
tities  are  <Rven  by  (C-100).  Table  C-7  gives  values  of  X  ,  ns  determined 

from  (C-165),  for  M'25"  and  y  =  I  . 

Table  C-7.  Thresholds  Ree.uired 
for  Processor  VII  for  M  25  and 


M  1;  Unknown 

Signal  Frequency 

IJF  (**  1) 

X(Pp,  M  25) 

io_1 

2.  1(1062 

_■> 

10  “ 

3.  2-1388 

io-:i 

•1.  10957 

10“‘ 

3.  0227 

Approximate  Detection  Probability 

The  RVs  u  and  V  in  (C-135)  are  approximately  Gaussian  for  large  M.and 
we  can  concentrate  on  their  second  moments.  We  have  already  evaluated  the 
mean  and  variance  of  u.  in  (C-121)  and  (C-121).  The  mean  of  RV  V  is  zero, 
as  may  be  seen  using  (C-ll(i).  The  mean  of  ItV  uv  can  be  determined  in  a 
manner  similar  to  that  in  (C-122)  and  (C-7  23);  the  result  is  found  to  bo  zero. 
The  mean-square  value  of  RV  V  is 

-  jAjt  ife  e[ (S„ c C, C.  5*.,)} 

+  J&C  [  Cs-r  f-«  5»2)(s„,  c„  -  CK.,S.| 

t  1e  {(%  Q., 

+  %-V5fe  T,  -  c  5,_$ | 


nT^rT-r~TrT?J;r^^  „r ..fM.T^.rT..r.v,  T.I-  ).tf,.i.jMr.«TTTi.KTi...  ■■» j.M^,uI«ii^iinaH«iMnii  u;ij.u«Mmii.x  uiffPiu  ink.  iwwiwrm<^ygfgra»BB^«*f»^^ 
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The  four  averages  in  (C-167)  are,  respectively, 

z??r  ?  € ,  -  S'  €,  o. 

Therefore  the  variance  of  RV  V  is 

The  joint  PDF  of  RVs  u  and  V  is  approximated  according  to 

f(“/v)  -  (2tf  r*) '  e*F  [‘  ^5^-  -  . 

The  probability  of  detection,  using  (C-134),  is 
Pp  ~  Pro\>(^  >1~)  ~  R^k(v/hV'7r  > 


iviv  (lVr„  rvY  ‘  exp[-  0  - 


where 


C-dxD-i]  F  fD/ 

Wi„ 

*;  ’ 


e  = 


T_ 


(M  -/)Z<T^'  f-y 


i-  ^  ■ 

uv 


(C-168) 


(c-:f‘9) 


(C-170) 


(C-171) 


(C-172) 


For  m=0  ,  changing  to  polar  coordinates  in  (C-171)  an  1  employing  (C-130)  and 
(C-1G9),  there  follows  for  the  approximate  false  alarm  probability' 

%  =  e^pl-xyz),  (C-173) 

where  normalized  threshold 


(C-174) 


92 
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which  agrees  with  (C-ii2)  and  ^C-156)  for  jX-o  ■  id  a  =  I  ,  respectively.  Then 
the  quantities  (C- 172)  that  are  necessary  tor  P  in  (C-171)  arc  expressible  in 
terms  of  X  as 


n  2r^  P(lt>^  V  ^  =  c=  ^d2,  cl=-^ 

a  77.  ^  ;  <V 


(C-175) 


The  quantities  mu  ,  <Tj^  ,  and  are  themselves  available  in  (C-121),  (C-12-1), 
and  (C-1G9),  respectively. 

Equations  (C-171),  (C-173),  and  (C-175)  constitute  the  approximations  to 
PD  and  Pr  for  arbitrary  jx  .  Threshold  X  is  determined  from  (C-173)  fora 
specified  and  substituted  in  (C-171)  to  evaluate  .  However,  for  la-0  or  jx~\  , 
X  should  be  determined  from  (C-141)  or  (C-1C>5),  respectively,  and  substituted 
in  (C-171)  for  evaluation  of  . 

Approximations  for  Email  m 

For  small  SNRs  per  sample,  vv>  «|.  The  approximations  in  (C-130)  are 
valid  in  this  case,  as  is 

^  [0)  /  P(l  +  /)  -  0-  .  (C-17(i) 


Then  (C-171)  becomes 


> 


duclv  (2-rrSyf)  exp 


(C-177) 


A5  hi  ->  0. 


-  «  y  rM  ^ 

where  we  <  hanged  to  polar  'ordinates  and  used  (C-8),  (C-17(i),  (C-171),  and 
(C-130).  The  factor  involving  p  in  (C-177)  was  encountered  previously  in 
(C-131)  and  was  tabulated  in  table  C-5;  the  comments  made  there  are  relevant 
here  also.  It  is  important  to  note  the  ha2 -dependence  in  (C-177)  versus  the 
linear  M -dependence  in  Processor  II,  (C-7),  and  in  Processor  V,  (C-G8).  Thus, 
in  addition  to  the  suppression  effect  caused  by  the  lack  of  knowledge  of  signal 
phas  j»0  ,  the  phase-difference  processor  suffers  further  small-signal  sup¬ 
pression,  for  unknown  signal  frequency,  for  all  values  of  weighting  constant  jK 
in  (C-133). 
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Appendix  D 

SIMULATION  PROCEDURE  AND  PROGRAM 
D.  I  SIMULATION 

In  order  to  determine  the  ROCs  tor  the  processors  in  tables  1  and  2,  10,000 
independent  trials  v/ere  ennducted  lor  each  processor.  A  trial  consisted  of  the 
generation  of  2M  independent  zero-mean  unit- variance  Gaussian  RVs.  To  keep 
computer  time  and  storage  at  reasonable  levels,  several  different  mean  values 
were  simultaneously  added  (in  different  summation  registers)  to  each  of  the 
RVs,  which  were  then  subjected  to  the  various  operations  shown  in  tables  1  and 
2.  The  number  of  threshold  excursions  of  the  resultant  RVs  were  evaluated  for 
a  wide  range  of  threshold  settings.  A  plot  of  the  relative  number  of  threshold 
crossings  for  a  particular  nonzero-mean  versus  the  relative  number  lor  a  zero 
mean  gave  one  curve  of  a  ROC.  A  sample  program  for  Processors  1-VI  is  given 
in  section  D.  2. 


For  a  given  added  mean  value,  the  number  of  excursions  versus  threshold 
settings  are  highly  correlated  RVs,  resulting  from  the  method  described  above 
for  generating  the  ROCs.  This  means  that  if  the  simulated  ROC  is  above  (below) 
the  true  ROC  at  one  value  of  ?F  ,  it  will  very  likely  be  above  (below)  the  true 
value  at  other  nearby  values  of  PF  .  The  way  we  have  chosen  to  control  this 
effect  is  to  take  many  trials,  namely,  10,000.  Thus  the  standard  deviation  of  a 
particular  probability  evaluation  is  (p(|-p)/io, ooq)^  s  ,  005"  ,  where  R  is  the 
probability  of  the  particular  event. 


In  order  to  determine  the  correlation  of  the  estimated  ROC  values,  consider 
the  following:  The  relative  number  of  excursions  above  threshold  A  for  /V  trials, 


where  U equals  I  if  ^>0,  mid  equals  0  if  0.  The  mean  of  RV  t(j\)  is 


rR)  =  UP- a)  -  I-  F(A)  *  QW, 


(R-2) 


where  ?(/)  is  the  probability  that  RV  x  remains  less  than  A  .  If  we  keep 
the  same  realization  of  the  RVs  as  in  (D-l),  but  subject  them  instead  to 

threshold  B  ,  yielding  RV 


r{$)  ^ 


<»-») 
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'  i~?r>  rnrv;,!^snfW'*fr",^">7  m 


**•  '.1'  «vu»|(w^ 


the  cross  correlation  coefficient  ot  >~(a)  and  r(®)  readily  follows  as 

Q  (wax  U,i$j)-  q(a)q(^ 

1°  [?(As)Q(a)?(-b)Q(b)]’/2 

using  the  independence  of  X„  and  XM  for  n^w.  So,  for  example,  if  yy-B,  (D-4) 
yields  p-  \  ,  as  it  must.  However,  if  A^O  and  B- I  ,  for  a  zero-mean  Gaussian 
RV  X„  ,  there  follows 

?(a)  ~  ,5\  P(s)  =  .Hl,  (D-5) 

1’hus  probability  estimates  as  different  as  .1300  and  .  841  are  still  rather  highly 
correlated.  In  order  to  avoid  overly  distorted  estimates  of  the  ROCs,  a  large 
number  of  independent  trials,  N  ,  must  be  used. 

D.  2  SAMPLE  PROGRAM  FOR  PROCESSORS  I 
THROUGH  VI 

Parameter  m=  21j-hi.-,s=iooo#  trials=ioooo,numsnr=s 

DIMENSION  ell :  UNS) ,P12(bINS) *P13(olNS) pP14(B1KS) .Plb(UlNS) 
P2l(iilNS)  *P22  (BINS)  *P23(ulNS)  pP24(DINS)  .P2b(BINS> 

1r  e31 (nINS) rP32(blNS) *P33(b!NS) pP34(B1NS> *P35(aINS) 

i*  e4  1  (til  NS)  *  P42  (BINS  )  *  P43(blNS)  »P44  (BINS)  *P45(  BINS) 

i.  Pbl(rilNS) *  P52 (BINS) *  P53 ( b I  NS ) *Pb4(BlNS) pPBb(BINS) 

i*  Pbl(BlNS) <  P62 ( BINS )  pP63(lInS) *P64(B1NS) *Pto5(8INS) 

DIMENSION  sA ( NUMSNk )  , So (NUMSNR) *  SC ( N( IMSNR ) * SD(NUMSnK) p  SF (NUMSNK) r 
iSF  (NUMSNR)  iA(o)  *  6  (  o  )  ,  Z (  200  )  *  R  (  b  )  *Se( NUMSNR) 

DIMENSION  el (Bl NS* NUMSNR) #02(BINS,NUMSNR) . 03 ( blNS * NUMSuR } , 

4  04  (  8  I  NS  i  ivUMSNR  )  f  05  (  H  I  NS  ,  NUMSNR  )  *  0b  ( B I  NS  *  NUMSNR  ) 

REAL  Ll  (NUMSNR)  *L2  (NUMSNR)  fL3(  NUMSNR)  >L4  (NUMSNR)  *  Lb  (  NUMSNk  )  * 

4Lo( NUMSNR) pMAXMpMAaMSQpMEANpMSTART 
EQUIVALENCE  (01 ( 1 p 1) .PI  1  ) 


EuUIVAlENCl 
EuUl  VAlENCl. 
EeUI VAlENCl 
EQUIVALENCE 
EQUIVALENCE. 
EuUl VAlENCl 
Eoul VALENCt 
EOUI  VALENCt. 
EOUl VAlENCl 
EuUIVAlEnCl 
equivalence 


(01  <  1  pH) 
(02(1*1) 
(02 ( 1 ph  ) 
( 03 ( 1 p  1  ) 
(  0  3  ( 1  f  4  ) 
(04 ( 1 p  1  ) 
(04 ( l pH) 
(OMIpI) 
( 05 ( 1 p  h ) 
(Oo(J.pI) 

(Opj  (  1  P  H  ) 


i  P  1  4  ) 
.P21  ) 
.P24) 
-  PA  l  ) 
'  P34  ) 
.PhD 
•  P44  ) 
.  Pb  1  ) 
.  P54  ) 
.Pbl  ) 


(01 ( 1 p2> ,Pl2 ) , (01 ( 1 p3) i P 1 3 ) 
(01 (1*5) ,P15) 

( 02  ( 1  *  2  )  p  P?  '; )  *  ( 02  (  1  *  3  )  pPt3) 
(02(1*5) *  P?5 ) 

(03(1*2) ,P32) * (03 ( 1 p3) *P33) 
(  0  3  ( 1  p  5  )  p  P  3  5 ) 

(04 ( 1 , 2) ,P4  2 ) p (04 ( 1 p  3) ,P43) 
( 04 ( 1,5) p  P 4 5 ) 

(05(1*?) pP52) r (05(1*3) pPa3) 
(05  (  US)  ,P55) 

( 06 • 1  .  2 ) , P62 ) 

(06(1*5)  P65) 


(06(1*3) ,Pb3) 


,Po4) 

OEFINE  F(I*K)  =  I*K+l(  1-SI6M1  *  I*K)  )/?)  *34359736367 

MSTART-.4 

DuLTAMz. 1 

K  -  b  *  *  1  S 

I=b28 1 

00M=1 , /M 

0uMSQ=00M**2 

OOT-D/TRIALS 


9(i 
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SuM-SQnT (OoM) 

MAXM-Dt.LT Am*  (imUMSNK-1  )  +MS1 ART 

MAXMSQrMAXN**2 

All) =  -4 . *SuM 

Bl 1 )zMAXM+4.*SQm 

0(2)  zMaXMSu+2  .  *00M  +  8,* •jQMfcbGR  T  (  MAXMc->0+0GM) 
A  (  J ) =  2  , -a .  »SOM 

B  (  j, )  z2  .  +MA  aMSO  +  R  .  *bQM*SQPT  IMAXM.S.Hl  .  ) 

A ( 4 ) z-2 • *SukT 12.) *SQM 
B(4):l, 

A  (  b  )  -  0  , 

B(b)  =  l , 

A(6)=A(3) 

Blb)zE3(3) 

Do  11  NPzl.b 
R (NP) =o!Nb/ (B(  .P)-A(NP)) 

DO  1  lT=lrlKlALS 
DO  9  Jr  1 » NuMSNK 

SA( J)=0. 

Sbl J)=0. 

SC  < J)=0. 

Sul J)=0. 

St(J)=0. 

SF (J)=0. 

SOI J)=0. 

CONTINUE 

PS I =6. 283 lo531 /TRIALS* IT 
CPrCOS  l  PS  I ) 

SP=SIN(PSI) 

DO  2  IM=1»h 
I rF ( I f  K ) 

W=T I NORM (Float ( I )/ 34 359738367. »$3) 

GO  TO  4 
PRINT  br 

format  i '  Problem* ) 

I=F(I*K) 

V=T  INOkM  (Fi_0AT  { I  ) /o435973b367 .  ,56) 

GO  TO  7 
PRINT  b » 

DO  8  Jr 1 , NuMSNR 

Me.  ANrDLLT  Am*  (  J-l )  +MS~'wT 

IF(J.Eu.l)  mEAi^-0 - 

WPMzto+ME AN 

XzwPM*CP-V*SP 

YzV*CP-s-WPM*SP 

SA  (  J  )  =SA  (  J  )  +  X 

Sb( J)=SB( J)+Y 

XS=X*X 

YS=Y*Y 

SC(J)=bC(J)+XS 
SUI J)=bD( J)+YS 
T  rbQRT ( XS+YS ) 

IMJ.Eu.l)  60  TO  41 

SG(J)=SG(J)+L0G(BSSL(T*MEAN,2)) 
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41  T=l./T 
bt-.(J)=bE(J)  +  X*T 
Sf-  <J)=i>F(J)+Y*I 

S  CONTINUE 

2  CONTINUE 

DO  10  J-1mnUMS'jP 
LI  (  J)=OOM*  (SA  (  J)  +LP+S0  (  J)  *bP) 
Lt(J)=00MSw*(SA(J)**2+bBJ J>**2) 
Lo(J)=00M*lSC(J)+^O(J) ) 

L4  (  J)  =OOM*  ( bE ( J )  +CM+SF  ( J ) +bP ) 

Lb ( J) =OOMSu*  (  Stl  <  J  )  **2+bF  (U>  ) 

IF (J.GT.l)  GO  TO  42 
L6< J)=l3< J) 

GO  TO  10 

42  MLAN=DELTAn*  ( J-3  )+>'.STAhT 
Lo( J)=4,*00M/MLAN**2*S&( J) 

iC  CONTINUE 

C  CAuCULATION  OF  PDF 
DO  12  u=l»NUMSiR 
N- ( L 1 ( j ) -A i 1 ) )*R(1) 

N-MAX  ( fj  i  1  ) 

N=NIN(N>ttlNS) 

01 (N* J)=01(N» Jl+OOl 
N=(L2(J)-A12) )*P(2) 

N=MAX(N»1) 

Nzm  I N  (  u  *  B  Ii-tb  ) 

02lNr J)=02lNr Jl+OOl 
N=(L3(J)-Al3)  )  +  P<  J) 

N=MAX (N> 1 > 

N=MIN(N»BInS) 

03(N, J)=03lNf Jl+OOl 
N=(L4( J)-A(4) ) *P(4) 

N-MAX (Nr  1 ) 

N=MIN ( N i BINS ) 

04 (N» J)=04(Nr J)+OUl 
N=(L5( J)-A(5) ) *R(b) 

N=MAX (N» 1 ) 

N-M IN  <  h  t  B I inS  ) 

Ob(N, J ) =OS  <  N  r J) +  0GT 
N=(L6(u)-Alb) )*'  (o) 

N=MAX (N» 1 ) 

N-M IN ( N  #  B 1 NS ) 

Ob(N. J)=06tN» J) +OUl 

12  continue 

1  CONTINUE 

C  CAl.CuLATIuN  OF  CDF 
DO  13  j-1 > NUMBnR 
DU  13  lB=2,BINb 

01  (BINb+l-lB-  J)=01  (BINl+1-IB,  J)  +01 l  < I Nb+ 2~ lb » J ) 
O^IblNu+l-lBf  J)=02(BINS+1-1B.  J 1 +02 ( .HNb+2-Ib > J ) 
0.1  (BINb+l“lB  »  J  >  =c  J  U3INb+l-lB»  J  )  +03  (olNb+2-Id  i  J  ) 

04  (dINbt-l-lB  f  J)  =04  (dINb+l-IBf  JI+04  li)lNb+2-Il'»u) 

0b(6INb+l-lB> J) =0b (BINb+l-lBf u > +0b ( 01 Nb+2-Ib » J) 
06  (  BINb+l-lB  »  J)=0o(BIN5+l~lBfJ)  +06  (BINb+2-IEi  i  J  ) 
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o  Continue 

CALL  MoUEbutZ.C) 

CaLL  SuB jto ( Z » 0 . , G . , 1 . , 1 , ) 

CALL  Ob JC To ( Z  » H 50. ,  1000  . *2850, ,2700 .  ) 
Call  LINESutZ, nlNS,Pll,P12) 

CALL  LlNt.SG(Z»i'lNS,Pll,P13) 
call  lineso  <  z  ..iIns.pii  ,  Pi4 ) 

CALL  LINESotZ.nINS.Pll ,P15) 
call  Subpli 

Call  LlNt,Sb(Z.olNb.P21*P22) 

Call  L I NEbu ( Z , B INs , P21 »  P23 ) 

Call  ulfvEbu  (  Z»  B I  NS , P21  »  P24  ) 

CALL  LINFSu ( Z , HINS.P21  ,P2b) 

CALL  buBPLl 

call  linesou*  'Ins.P31 ,P32> 

Call  L  I  NESu  (  Z  » •'  I  Nb ,  P31 ,  P33  ) 

CALL  LINESutZ, . ilNS, P31 ,P34) 

CALL  LI  NESutZ.nl NS »P31  ,  P3b ) 

Call  SuBPLI 

C«Ll  LINESutZ, <iINS.P4l,P42) 

Call  LINESutZ,  1 N  b  *  P 4 1 .  P  4  u  ) 

Call  LINESutZ, ;.tlNS,P41,P44) 

CALL  LINESutZ.. 1INS.P41  ,P4b) 

CALL  bLiBPLl 

CALL  L iNESo ( Z  ,  u I Nb , Pbl , Pb2 ) 

Call  LINESutZ,  'lNb,Pbl.Pb3) 

CALL  LINESutZ.  B I  NS , P51 , P54 ) 

CALL  UNtSotZ.  >INS»P51  ,P55) 

CALL  bUMPLI 

CALL  LlNLbO(Z,!tINb.P61  ,  P62  ) 

CALL  LINESutZ, MlNS.P6l.P63) 

CALL  LlNEbutZ,HlNS,P6l ,P64) 

CALL  LINE  j(Z,M INS, P61.P65) 

CALL  SuF' 

PE  L  = ,  1 

DO  21  iBzl.BINS 
KUIB 

c-l  IE  (Pll  (  IB)  .LT.PPL)  GO  TO  22 
t2  DO  23  1B=1  .BIN-5 

K  ^  -  I B 

c3  IE (P21 ( IB) .LT.PF L)  GO  TO  24 
c. 4  DO  25  IB=1  .BIN. 

K3=1B 

--o  IE  (P31  (  IB)  .LT.PPL)  GO  TO  26 
c. o  Do  27  IB=  1  .BIN  > 

K4zIB 

<-!  IE  (P41  t  IB)  .LT  .PEL)  GO  TO  26 
£6  DO  29  i  1 , B I Nb 
Kb=IB 

IMP51  (  IB)  .LT.pFL)  GO  TO  30 
3c  DO  31  1B=1  ,  Ee  I  Nb 
Kt-  I B 

ui  IE (Pbl ( IB) .LT.PEL)  GO  TO  32 
32  CALL  SuUJEOtZ.0.  ,0.  ,-pFL,  1 .  ) 

CALL  LlNESo (Z.BINS+1-K1  .PlltKl)  ,P12(Kl) ) 
Call  LlNESoiZ,"lNb+l-Kl , PlltKl), P13(K1)) 
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CALi.  LlNEbo(Z.l« 

CALL  L  INLSli  (  Z  .  1  *  .  1  »  0  ,  ) 

CALL  L 1 NESo ( 2 . 1 ' 0 >  .  0  ,  ) 

CALL  L INESu ( Z » 1 » • 1  * • 1 ) 

CALL  SLTSMu(Z>30,  ,‘j) 

Du  15  lLzl.9 

CAlL  L1NESo(Z»C. .01.il. 0.) 
CALL  LINE,  bo  (Z»1*.0'*1L*1<) 
CALL  LINLSo <Z . 0 . 0 . # . 1*1L> 
CALL  LINESo ( L »  1 »  .  1  *  •  1*IL) 
CALL  PAGEGtZ.CI.2fl) 

CALL  SLTSMo l  Z  »  A  0  *  1  .  ) 
RLTURN 

end 
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